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MULTIVARIABLE g-RACAH POLYNOMIALS 



J. F. VAN DIEJEN AND J. V. STOKMAN 



Abstract. The Koornwinder-Macdonald multivariable generalization of the 
Askey- Wilson polynomials is studied for parameters satisfying a truncation 
condition such that the orthogonality measure becomes discrete with support 
on a finite grid. For this parameter regime the polynomials may be seen as 
a multivariable counterpart of the (one-variable) q-Racah polynomials. We 
present the discrete orthogonality measure, expressions for the normalization 
constants converting the polynomials into an orthonormal system (in terms of 
the normalization constant for the unit polynomial), and we discuss the limit 
q — » 1 leading to multivariable Racah type polynomials. Of special interest 
is the situation that q lies on the unit circle; in that case it is found that 
there exists a natural parameter domain for which the discrete orthogonality 
measure (which is complex in general) becomes real-valued and positive. We 
investigate the properties of a finite-dimensional discrete integral transform 
for functions over the grid, whose kernel is determined by the multivariable 
(jr-Racah polynomials with parameters in this positivity domain. 



1. Introduction 



Some y ears ago, Koornwinder extended a construction of Macdonald |M1 



(see also | VK] |) to arrive at a multivariable generalization of a family of basic 
hypergeometric polynomials commonly known as the Askey- Wilson polynomials 
| AW 2 , GR | . The multivariable polynomials of interest depend rationally on a num- 



ber of parameters and for parameter values in a certain domain they form an 
orthogonal system with respect to an explicitly given (positive) continuous weight 
function with support on a (real) n-dimensional torus (where n denotes the num- 
ber of variables) . Recently, it was shown that the parameter domain for which the 
multivariable Askey- Wilson polynomials admit such an interpretation as orthog- 
onal polynomials may be extended if one allows the corresponding orthogonality 
measure to have a partly continuous and partly discrete support |^j. (Thus further 
generalizing the corresponding situation in the case of one single variable, where 
the phenomenon of discrete masses emerging in the Askey- Wilson orthogonality 
measure was already known to occur [ AW2[ .) 



In the present paper we will demonstrate that for a different parameter regime 
satisfying a certain truncation condition, the multivariable Askey- Wilson polyno- 
mials can be reduced to a finite-dimensional orthogonal system with respect to a 
purely discrete weight function living on (i.e. supported on) a finite grid. The 
polynomials amount for these parameters to a multivariable generalization of the 



q- Racah polynomials [ AW1 , GR |. We will compute the normalization constants 
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turning the polynomials into an orthonormal system with respect to the discrete 
orthogonality measure (in terms of the corresponding normalization constant for 
the unit polynomial), and also discuss the limit q — > 1 giving rise to multivariable 
Racah type polynomials. 

In general we will not worry much about the positivity of the weight function 
and we will, in fact, for most of the time allow parameters to be such that the 
discrete orthogonality measure becomes complex. However, as will be outlined 
below in further detail, it is possible to restrict the parameters for \q\ = 1 in a 
rather natural way to a subdomain such that the discrete orthogonality measure 
for the multivariable g-Racah polynomials becomes a genuine positive measure. For 
parameters in this positivity domain the restriction of the multivariable g-Racah 
polynomials to the grid points entails (after renormalizing) an orthonormal basis 
for the finite-dimensional Hilbert space consisting of all (complex) functions over 
the grid. From a functional-analytic point of view, the renormalized polynomials 
determine the kernel of a unitary finite-dimensional integral transformation for 
these grid functions. 

The material is structured as follows. In Section || we first recall the definition 
of the multivariable Askey- Wilson polynomials for generic parameters as eigen- 
functions of the Koornwinder-Macdonald second order analytic ^-difference oper- 
ator. Next, the discrete orthogonality measure for the polynomials is introduced 
in Section |^ and it is explained how the multivariable Askey- Wilson polynomials 
descend to a q-Racah type finite-dimensional orthogonal system when being re- 
stricted to the finite grid. Crucial in the orthogonality proof is the observation 
that the Koornwinder-Macdonald second order g-difference operator diagonalized 
by the polynomials is symmetric with respect to the discrete inner product (just as 
this observation turned out to be essential when dealing with a purely continuous 
or mixed continuous/discrete orthogonality measure jM| g §). The proof for the 
symmetry property of the g-difference operator (in the discrete context) is relegated 
to an appendix at the end of the paper (Appendix |A]) . The orthonormalization con- 
stants are given in Section ^. Their computation, which is outlined in Appendix |b| 
makes use of a recently introduced system of recurrence relations (or Pieri type 
formulas) for the multivariable Askey- Wilson polynomials [D2, D3]. In Section || 
we continue by discussing the transition from the basic hypergeometric level to 
the hypergeometric level (q — > 1). In this limit — which one might also interpret 
as a transition from trigonometric polynomials to rational counterparts obtained 
by sending the period of the trigonometric functions to infinity — our multivariable 
q- Racah type polynomials degenerate into multivariable Racah type polynomials. 
The paper is concluded in Section ^ with the characterization of a parameter domain 
for which the weights determining the orthogonality measure for the multivariable 
g-Racah polynomials become positive when q lies on the unit circle, and a descrip- 
tion of some properties of the finite-dimensional discrete integral transform for grid 
functions that is associated to the polynomials with parameters in this positivity 
domain. 



Note. Most objects of interest in this paper (such as the polynomials, weight func- 
tions, difference equations, normalization formulas, etc.) depend rationally on a 
number of parameters. Below it will always be assumed (unless explicitly stated 
otherwise) that the parameters are such that one stays away from singularities (this 
is the case generically) , without repeatedly stressing this point each time. 
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We will be mainly concerned with the study of complex orthogonality properties 
of our polynomials. We shall say that a basis {pi}i for a complex vector space H 
endowed with a nondegenerate bilinear form (•, •} is orthogonal (with respect to 
(•, ■}) if (pi,pj) — for i 7^ j. Furthermore, an orthogonal basis of TC will be called 
orthonormal (with respect to (■,■)) if in addition the normalization is such that the 
quantities {pi,Pi} are all equal to one. 



2. MULTIVARIABLE ASKEY-WlLSON POLYNOMIALS 

In this section the definition of the multivariable Askey- Wilson polynomials as 
eigenfunctions of the second order Koornwinder-Macdonald q-difference operator 
is recalled. The approach followed here is very much in the spirit of Macdonald's 
original treatment in Ml]. 

Let 7i = C[zi, . . . , z n , z^ , . . . , z" 1 ] be the space of Laurent polynomials in the 
variables zt, . . . ,z n . On this space the (BC type) Weyl group W = S n K (Z2)™ 
(i.e. the semidirect product of the permutation group S n and the n-fold product 
of the cyclic group Z 2 ) acts naturally by permutation and inversion (zj — * zj 1 ) 
of the variables Z\,... ,z n . The subspace TC W of Ti. consisting of the VF-invariant 
(Laurent) polynomials is spanned by the symmetrized monomials 

(2.1) m x {z)= AeA ' 

where A(= W \ Z n ) denotes the integral cone 

(2.2) A = {A e Z" I Ai > A 2 > • • • > X n > 0} 



and the summation in (2.1) is meant over the orbit of A under the action of the 
group W (which acts on vectors A S Z" by permuting and flipping the signs of 
the vector components Ai, . . . , A„). For future refer enc e we will partially order the 
basis elements m\, A € A, by defining for fj,, X £ A (2.2) 

(2.3) fi<X iff A J for k = l i---i n 

l<j<k l<j<k 

(and n < A if fi < X with ^ X). 

In |H , Koornwinder introduced the following generalization of the second order 
(BC type) Macdonald q-difference operator 



(2-4) D= J2 faz) (T hq 1) + V.j(z) (Tr 1 - 1)), 

with 



l<j<n 



(1 - to25)(l - ti4)(l - t 3 2f)(l ~ hz 



(l-zf){l-qzf) 

(\-tz]z k )(\~tz]z- k ^ 

11 (1-^^)(1-^^) ' 
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and showed that this operator is triangular with respect to the partially ordered 
basis of monomial symmetric functions: 



(2.5) 



E ^ 



with £7 X eC[5 ±1 ,i,t 0) ti,i 2 ,i3] 



(i.e. the expansion coefficients (or matrix elements) E\^ depend polynomially on 
t and to, . . . ,£3). The leading coefficient (or diagonal matrix element) -Ea.a in 
(fl|) reads explicitly §§ 



(2.6) 



£a,a = 



E 

l<j<n 



9 1 t tit 2 t 3 t 



2n-i-l /„A 



-l/„-A, 



-1) + ^(g 



!))• 



The triangularity (2.5) of 13 implies that the eigenvalue problem for the g-difference 
operator in the space TC W is essentially finite-dimensional, because it can be re duc ed 
to the invariant subspaces of the form H.Y = SpanjTO^J-^gA^A (with A 6 A Q2.2|)) . 
The eigenvalues of D in Tf^ are given by the diagonal m atri x elements E\\ ( |2.6[ ), 
A 6 A. It is immediate from the explicit expression in (|2.6| ) that E^ ^ E\\ as 
a polynomial in the parameters q^ 1 , t and to, ■ ■ ■ , £3 if \i ^ A. In other words, the 
eigenvalues are nondegenerate as (Laurent) polynomial functions of the parameters. 
The Koornwinder-Macdonald multivariable Askey- Wilson polynomials are now — by 
definition — the corresponding eigenfunctions p\, A 6 A. 

Definition. The multivariable Askey-Wilson polynomial associated with a (domi- 
nant weight) vector A £ A {2/2) is the (unique) monic VF-invariant Laurent poly- 
nomial of the form 



(2.7a) p x (z)=m x (z) + 



such that 



E 



cx.^m^z) with £ C(q, t, t , ti, t 2 , h), 



(2.7b) 



Dp\ = E X \P\- 



Following Macdonald [MI] (see also [D2, D3 



it is possible to write down a 
somewhat more constructive formula for the polynomial p\ in terms of the mono- 
mial m\, the operator D and the eigenvalues E^ (// < A) 



(2.8) 



Px 



n 



D-E„ 



Ex, 



E„ 



(Notice that the r.h.s. is well-defined as a rational expression in the parameters in 
view of the fact that the denominators are nonzero as (Laurent) polynomials in the 
parameters.) The validity of this representation fo r th e multivariable Askey-Wilson 
polynomials is easily verified by inferring that p\ ( |2.8| ) satisfies the defining proper- 
ties ( [2. 7a ) and ( 2.7b ). That the r.h.s. of (2.8) is of the form in ( 2.7a ) is clear from 



the triangularity of D; that it also satisfies ( 2.7b ) is immediate from the fact that the 
operator n^eA ^i<a(^ — annihilates the subspace TiY — Span{m^}^eA,/i<A 

as consequence of the Cayley-Hamilton theorem (and hence applying (D — E\,\) 
to the r.h.s. of (2.8) yields zero, which is precisely (2.7b)). 

In the one variable case, there is no i-dependence and the eigenvalue equation 
Dp\ — E\_\p\ amounts in that case to the second order g-difference equation for 
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the Askey- Wilson polynomials. Thus, the polynomials p\ then reduce to (monic) 
Askey- Wilson polynomials |AW2| , [GR 

, n n x , (toti,tot2,tot3;q)x fq~ x , hhhhq toz, tez^ 1 

(2.9) px{z) = A , , . . , , ;q,q 

toKhhhhq ,Q)x V toti, tote, tote 



with A = 0, 1,2, . . . Here we have employed standard notation (see e.g. AW2 , GR|) 
for the basic hypergeometric series 

, ( ai,...,a s+1 \ _ (ai,... ,a s+ i;q) m z m 
S+10S l ^---A ,q,Z )~£- Q (h,... ,b 8 ;q) m (q;q) m 

and the g-shifted factorials 

(oi,... ,a s ;q) m = (ax;q) m ■ ■ ■ (a s ;q) m , 

(a;q) m = (1 - o)(l - aq) ■ ■ • (1 - aq m - 1 ) 

(with (a;q) = 1). 

3. A Discrete orthogonality measure 

It follows from jKj] that if the parameters satisfy the constraints 

(3.1) 0<g<l, -Ki<l, |i r | < 1, r = 0,1,2,3, 

with possible non-real parameters t r occurring in complex conjugate pairs and 
pairwise products of the t r not lying in the interval [l,oo), then the polynomials 
px, A G A constitute an orthogonal system with respect to a continuous weight 
function A AW : 



(3.2) /•••/ px(e t3: )p fl (e i *)A AW (x)dxi---dx n = Q if A ^ ,i 

where 



(3.3) A AW (a;) = J[ 



^ e i{e l x ] +e 2 x k ). q ^ oo 



11 (ie i ( ei ^+ e2a:fc );g) 00 

£1 ,E2=±1 

x -q (e 2ie ^;g)oo 



(to e tex i , t\ e iex i , te e %eXi , te e tex i ; q) c 

1 < j < n 

e=±l 



Here we have used the notation e lx = (e lXl , . . . , e lXn ) and (cf. above) 

oo 

(ai, ... ,a r ; q)^ = (a^q)^ ■ ■ ■ (a r ; q) x , (a; q)^ = JJ (1 - aq m ). 

m=0 

It should be noted that at this point of our presentation it is not even obvious that 
the multivariable Askey- Wilson polynomials defined in the previous section actually 
exist for all parameter values in the above domain, since the domain might a priori 
contain parameter values for which the coefficients ca,^ (2.7a) have a singularity. 
However, that such singularities indeed do not occur is seen from an alternative 
characterization of the multivariable Askey- Wilson polynomial px as the polynomial 
of the form 

(3.4a) px(z) = mx{z) + 2J c A:AI m M (z) with c A , M e C, 

peA, fj,<X 
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satisfying 

(3.4b) /••■/ Px (e lx )m fl (e™)A AW (x)dx 1 ---dx n = for /j < A. 



— "7T J —TT 



(In [|k| these two properties were in fact used to define the multivariable Askey- 
Wilson polynomials.) It is clear that the polynomials determined by the conditions 
( [3.4a ) , ( 3.4b| ) are well-defined for all parameter values in the domain given by ( |3.l[ ) 



and, furthermore, that they are continuous in t and to, ■ ■ ■ , £3 for these parameter 
values. By showing (as was done in 0) that the difference operator D ( |2.4[ ) is 
symmetric with respect to the L 2 inner product with weight function A AW , i.e., 
that 



— 7T J — 7T 



{Dm x )(e iX ) m M (e ia: )A AW (x)dxi ...dx n 



m x {e lx ) {Dm f ,)(e ix )A AW (x)dx 1 ■ ■ ■ dx n 
and combining this with the triangularity property in (|2.5|), one finds that the 



polynomial p\ (z) of the form (3.4a), (3.4b) satisfies the eigenvalue equation Dp\ 



E\,\P\, which shows that for the parameter domain of interest this alternative 
characterization leads us to the same polynomial p\ as in the previous section 
(and hence that the coefficients are regular for these parameter values). The 
orthogonality ( |3.2| ) of the multivariable Askey- Wilson polynomials p\, A £ A with 
respect to the weight function A AW ( |3.3| ) now follows, first for generic parameters in 
the above domain from the fact that they are eigenfunctions of a symmetric operator 
D corresponding to different (real) eigenvalues, and then for all parameters in this 
domain by a continuity argument jKJ]. 

Very recently, it was shown that the parameter domain for which the multivari- 
able Askey- Wilson polynomials admit an interpretation as orthogonal polynomials 
may be further extended. Specifically, it follows from ^ that for t — q m , with 
to an arbitrary nonnegative integer, one can remove the constraints \t r \ < 1 in 



(3.1) (while still keeping the other restrictions on these parameters though) to end 



up with an orthogonality relation for the polynomials p\, A £ A consisting of a 



continuous part of the form (3.2) and additional mixed continuous, mixed discrete 
parts. 

The main purpose of the present section is to demonstrate that for parameters 
satisfying a truncation condition the multivariable Askey- Wilson polynomials give 
rise to a finite-dimensional orthogonal system with a purely discrete weight function 
living on a finite grid. 

First it is needed to introduce some notation. Let N be an arbitrary nonnegative 
integer and let An be the alcove of dominant weight vectors of the form 

(3.5) A N = {A e A | A < Nut} with uj = e% H h e„. 

(Here and below ej represents the jth unit vector in the standard basis of K™.) We 
will show that for parameters subject to a truncation condition of the type 

(3.6) tatbt"- 1 = q- N , 

with a,b € {0,1,2,3} (fixed) and b ^ a, the multivariable Askey- Wilson polyno- 
mials p\, A G An form a finite-dimensional orthogonal system with respect to a 
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(generally complex) discrete measure with support on the grid points 

(3.7) rq" = [r iq v \. . . , r^V . . , r n q^), v e A N 
where 

(3.8) Tj =t n ~ j t a , i = l,...,n. 

Specifically, the (complex) orthogonality relation becomes for these parameters (i.e. 
satisfying condition (|3.(j)) 

(3.9) Yl Px{rq u )P^)^ K {v) = for (A, p E A N ), 
i^eAjv 

where the weights A qR (i/), v £ Ajy are given explicitly by 

(3 - 10) AqR( ^ = c W Wr 



with 



(3.11) CI {v) = c (v) M — -2 

1 1 ( <t j r fc ; 9)^ + v» (*7j r fe ; q) Vj _ ^ 



l<.j<fe<n 

>< n 

l<j<n 

(3.12) C«» = coM [] ( 

l<j<fc<n 



{Tjiqhvj 



X 



n 



^ilo<r<3(' rT 3'9)i'3 / 

(qTjT k :q)v j+ „ k (qTjT^ 1 ^)^—, 
'llo<r<3(*r Vi;«0r/A 



1<J<„\ te^K- 



and 

(3.13) co(i/)= [] (^(ioiit 2 t 3 ^ 1 ) 1/2 )^ • 

l<j<n 

Here we have introduced discretized functions C±*{v) that are reminiscent (upon 
dualization) of the c-functions of Harish-Chandra in harmonic analysis (see e.g. 



He, H£]). 



Lemma 3.1. a. For v G A (|J) the functions C q ±{v) ( |3.1l| ), ( |3.12| ) are well- 
defined and nonzero as meromorphic expressions in the parameters t, to,... ,t% 
and q. 

b. For v e An (H)) the functions Cf{v) < ^u\) , ( ^12| ) are well-defined and 
nonzero as meromorphic expressions in the parameters t, to, . . . ,t% and q subject 
to the truncation condition ( [3.6| ). 

Here and below when stating that an expression is meromorphic or rational in the 
parameters i, to, . . . , £3 and q subject to the truncation condition ( ft.6| ), it is meant 
that after elimination of 4 (or t a ) with the aid of the relation t a tb = t 1 ~ n q~ N , the 
resulting expression is meromorphic/rational in the remaining parameters t a (or 
tb), tc,td,t and q (where t c and td denote the two parameters complementing t a 
and t b in {to,*i,*2,*3} such that {t a , t b , t Cl t d } = {t a ,ti,t 2 ,t 3 }). 

The proof of this lemma is immediate from inspection of the explicit expressions 
for CJ R (^) given above (no factor in the numerator or the denominator of Cj^'(i') 
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becomes identical to zero). Observe also that the functions CJ R (^) are rational in 
the parameters except for the square roots appearing in the common factor cq(i>) 
( 3.13| ). In the discretized weight function A qR (^) ( |3.1C ) these square roots can 



be collected (there is an even number of them) and rationality becomes restored. 



Notice also that for parameters subject to the truncation condition C} {v) (3.11) 
would be infinite if v G A \ A^r (i.e., when v lies in the cone A ( |2.2| ) but outside 
the alcove A^r (|3.5[)). In that case we have that v\ > N and hence the factor 



(tfeTi; q) Vl — {t a t b t n 1 ; q) vi in the denominator of C qR (^) (|3.11 ) is identical to zero 



when the truncation condition (3.6) holds. As a result, the weight function A qR (z/) 
( |3.10| ) vanishes on A \ Ajy for parameters subject to the truncation condition. 



The orthogonality relation (3.9) should be read as a set of identities for the 



multivariable Askey- Wilson polynomials p\{z), A G An with parameters subject to 



the truncation condition (3.6). It is clear from Macdonald's representation in (£ 
and the explicit formula for the eigenvalues E\.\ ( |2.6| ) that the polynomials p\{z) 
with this restriction on the parameters are well-defined as rational expressions in the 
parameters subject to the truncation condition (no denominator in ( |2.8| ) becomes 
identical to zero after imposing the truncation condition (|3.6|)). Since the same is 



true for the weights A qR (^) ( 3.1 0| ) , the terms in the l.h.s. of ( |3.9| ) make sense as 



rational expressions in the parameters q, t and to,ti,t 2 , t% subject to condition (3 

In the case of one single variable (n = 1) the orthogonality relation (^^) reduces 
to the well-known orthogonality relation 



E 

0<v<N 



(3.14) X P\W)p»M)&*{y) = 

for (A, M G{0,... ,N}), 



with 

(3.15) A qI » = 

for the (monic) Askey- Wilson polynomials given by p \(z) ( |2.9|) with parameters 
subject to a condition of the type t a tb = q~ N (b ^ a) [AW1, pR[ . In the discrete 



(1 — t 2 a c t 2u ) (tptqi tit a , t2t a ,t3t a ; q) v 

(tohht-jq- 1 )^! - t 2 a ) {t^qta, t^qta, *a Va, Va5 q) t 



case, with parameters subject to the truncation condition, these one-variable poly- 
nomials are usually referred to in the literature as q-Racah polynomials rather than 
Askey- Wilson polynomials. 

In the above formulas (and throughout the paper) t a and tf, denote two fixed, 
distinct, but otherwise arbitrary parameters from the set {to, t\, t$\. Since the 
monic (multivariable) Askey- Wilson polynomials p\(z) are symmetric with respect 
to permutations of the parameters to, . . . ,t% (this is immediate from their definition 



as monic eigenfunctions of the difference operator D (2.4)), the actual choices of t a 



and tb are all equivalent up to permutation of the parameters to,ti,t2,t^. Notice, 



however, that the choice of t a does influence the position of the grid points rq 1 " (3.7). 
In the case of one single variable one usually works with a different normalization 
of the g-Racah polynomials that breaks the permutation symmetry and designates 
a preferred role to the parameter to- (Instead of employing the monic polynomials 
of ( |2.9| ) one then omits the constant factor in front of the terminating 4^3 series.) 
It is in that case custom to take t a to be to (thus fixing the grid points to toq u , 
v = 0, . . . , AT). This then leaves three possible choices for the other parameter tb 
entering the truncation condition t a tb = q~ N , which leads one to distinguish the 
three cases toti = q~ N , t^ = q~ N and to^3 = q~ N ■ Here instead of treating each 
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of these cases separately we prefer not to make such distinctions and have tried 
rather to emphasize the symmetry in the parameters to, . . . ,t^ with our notation. 



In order to prove the discrete orthogonality relations (3.9) and analyze their 



properties in more detail we will need some further notation. Specifically, let us 
introduce a bilinear form (•, on the finite-dimensional subspace (of Ti. w ) 

(3.16) = Span{m A }AeA J v 
determined by 

(3.17) (f,9)f= J2 f(rql9(rqn^ R W) 

i^eAjy 

The orthogonality of the multivariable Askey- Wils on (or g-Racah) polynomials with 
respect to the discrete weight function A qR ( |3 .10 ) can now be phrased in terms of 
the following theorem. 

Theorem 3.2. For parameters q, t and to, ■ ■ ■ ,t% subject to the truncation condi- 
tion 

tahf 1 - 1 =q- N b^a 

( where N denotes an arbitrary nonnegative integer and a, b are two fixed, distinct 
but otherwise arbitrary numbers from the set {0, 1,2,3},), the multivariable Askey- 
Wilson polynomials p\(z), A £ An (|3.5|) satisfy the orthogonality relation 



(3.18) (PA,rf = for X + v (X,fieA 



N 



with the bilinear form (•, being defined by ( 3.17| ). 



The proof of the orthogonality theorem operates along the same lines as the 
orthogonality proof in Q for the continuous case (cf. the beginning of this sec- 
tion) and is based on the statement that (for parameters satisfying the truncation 
condition) the difference operator D ( |2.4| ) is symmetric with respect to the bilinear 
form (•, ■) N X . 

Proposition 3.3. For parameters in accordance with the assumptions in Theo- 



form (■, -) q N ( |3-17| ), Le - 



rem \3. 4 the difference operator D (2.4) is symmetric with respect to the bilinear 

/ N 



(3.19) (Df,g)f = (f,Dg)f (f,g eUf). 



The proof of this proposition has, in turn, been relegated to Appendix ^ at the 



end of the paper. The orthogonality (3.1q) now follows from the fact that the poly- 



nomials p\, A £ An are eigenfunctions of a symmetric operator D corresponding 



to eigenvalues E\ \ (2.6), A £ Ajv that are nondegenerate (i.e. distinct for different 
A) as rational functions in the parameters q, t and to, ■ ■ ■ ,t^ subject to condition 
&■ 

An important property of the bilinear form (•, is that for generic parameters 
(subject to the truncation condition) it is nondegenerate on the space W^f (3.16) 
(i.e., if for a certain / £ Tiff one has that (/, g)ff = for all g £ TCf^ , then / must 
be zero). 



Proposition 3.4. The bilinear form (•, -) q N (3.17) is nondegenerate on the space 



Ti. q j^ ( 3.16| ) for generic parameters q, t and to, ■ ■ ■ ,t% (subject to condition ( |3.6| ) ). 
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The proof of Proposition 3.4 readily follows after recalling that the weights 
A qR (^), v £ An are nonzero as rational expressions in the parameters subject 
to the truncation condition (no weight becomes identical to zero) and combining 
this with the following lemma. 



Lemma 3.5. For generic q, t and t a every function in the space TL q ^ (3.16) is 
uniquely characterized by its values on the grid points rq x , X E Ap 



IN. 



To prove Lemma 3.5 one uses that the determinant of the matrix 



(3.20) 



M 



(where the columns and rows ar e or dered by means of some (any) tota l order of the 
weight vectors in the cone A (2.2) extending the partial order (2.3)) is nonzero 
as a Laurent polynomial in q for generic t and t a (t his follows from the domi- 
nant behavior of the determinant for q — > +oo, see [D2, Lemma 4]). Hence, if 



J^ueAjv c ^ m ^( T< l v ) = for all v € An, then the coefficients must all be zero 



Since the dimension of the space of all functions over the grid (3.7) is equal to 
the number of points in the grid, which is for generic q, t and t a equal to the number 
of points in the alcove An ( |3.5| ) and hence identical to the dimension of the space 
(3.16), Lemma 3.5 actually says — when read in an appropriate way — that for 
generic q, t and t a any complex function defined on the grid points rq x , X G A^ may 
be identified with the restriction of a function in the space Tiff to the grid. In other 
words, for such generic q, t and t a every function on the grid can be approximated 
exactly by (the restriction to the grid of) a unique function in . 

Proposition 3.6. For generic q, t and t a any complex junction over the grid with 
points Tq v , v £ An can be represented exactly by the restriction to the grid points 



of a unique function in the space "hCS ( 3.16| ). 



In principle there is no real reason why the definition of the bilinear form 



'N 



(3.17) should be confined to the subspace TC^ and could not, e.g., be extended 
to the whole space TL W (which consist of all W-invariant Laurent polynomials in 



n variables). In fact, this is why the condition f.g G Tiff in (3.17) has been put 



between parenthesis. It turns out that Theorem |3.2| and Proposition 3.2 remain 
valid (with the given proofs applying verbatim) if one allows A and /x to be any 
dominant weight vector from the cone A (2.2) and / and g to be any function in the 
space Tt w , respectively. However, since such a generalization of Theorem 3.2 would 
imply that p\ is o rthogonal to TCj^ if A is not in Ajv, it follows from combination 
of Proposition 3.4 and Proposition 3.6 that in that case p\ must actually be zero 
on the grid. 



Proposition 3.7. For parameters subject to the truncation condition (3.6) one has 
that 



(3.21) 

ifX G A\ A N . 



V\{rq v ) = for v G A 



A* 



The upshot of Proposition 3.7 is that the multivariable Askey- Wilson polynomi- 
als p\(z), A G A with parameters subject to the truncation condition ( |3.6[ ) descend 
to a finite-dimensional orthogonal system spanned by p\(z), X G An when being 
restricted to the grid points rq", v G An- Furthermore, the results of this section 
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can now be alternatively summarized in the following way: for generic q and pa- 
rameters subject to the truncation condition, the restriction of the multivariable 
Askey- Wilson polynomials p\{z), A G Ajy to the grid points rq" , v e An yields 
an orthogonal basis for the (finite-dimensional) space of functions over this grid 
endowed with the bilinear form X^eAjv f( T 1 l ')9( T( l^)^ qK ( 1 ')- 

Remark 3.1. It is convenient to think of the discrete weight function determined 



by the weights A qR (^) (3.10) as arising from a continuous weight function A qR (z) 
restricted to the grid points rq v ( |3.7D (such that A^^rq") — A qR (^)). In principle 
such a function A qR (z) of course exists and is a priori highly non-unique. For 
< \q\ < 1 a possibility would be to take 

A qR (z) = V Y[ (t^^Hohhtsq- 1 )- 10 ^/ 10 ^ 
l<j<n 

(3.22) x n (d - WW - Z^) ^M^w?, q )oo \ 

l<j<k<n\ (tZjZk,tZjZ k ,q)co J 

x tt ( / 1 .. 2 1 <i z j 1 ggj -^2 1 i z j 1 izj ; q) oo \ 

i<j<rA 3 (t Zj,tiZj,t 2 z j ,t 3 Zj;q) oo /' 

where T>o denotes a normalization constant determined by the requirement that 



evaluation of A qR (z) in z = t = (n, . . . , Tj, . . . ,t„) — i.e., z = rq v (3.7) with 
v = — should yield the value one (since A qR (0) = 1). 

Remark 3.2. It is instructive to view the grid points ( |3.7|) on which the discrete 
masses of the orthogonality measure determined by A qR are positioned as being of 
the form q p+v , v G An, where p denotes an n-dimensional vector that is related to 
r ( |3.8| ) by a logarithm 

(3.23) p= (pi,... ,p n ) with pj =log(Tj)/log(g), j = l,...,n. 

Thus, up to an exponential scaling (z — > q z ) the grid consists of the points in the 



alcove An (3.5) translated over the vector p. 

4. Orthonormalization 

In this section we will present the normalization constants turning the multivari- 
able Askey- Wilson polynomials p\(z), A £ A^r (|3.5|) — with parameters subject to 



the truncation condition ( |3.6| ) — into an orthonormal basis for the space T~t N (3.16) 

vqF 
>N 



with respect to the bilinear form (•, (3.17). To this end it is needed to compute 



the quantities (px,px)'^, A e Ajv- The evaluation of these sums pivots on a pre- 
viously introduced system of Pieri type recurrence formulas for the multivariable 



Askey- Wilson polynomials |D2, D3|. At the time of their introduction the same for- 



mulas were also used to verify the explicit expressions (conjectured by Macdonald) 



for the squared norm of P\(z) with respect to Koornwinder's inner product (3.2). 
In both cases (i.e. the discrete and the continuous case) the mechanism leading 
to the solution of the orthonormalization problem is very similar. Hence, wc will 
refrain from providing full details here and refer the reader to Appendix |b| where 
the main ingredients of the recipe are outlined with an emphasis on those points at 
which the discrete case differs from the continuous case. 

We would like to add that exactly the same approach may also serve to obtain 
the orthonormalization constants (in terms of the norm of the unit polynomial) for 
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the parameter regime with the mixed continuous/discrete orthogonality measure 
considered in ^ . The answer will in that case be formally identical to that in the 
case of Koornwinder's purely continuous orthogonality measure A AW ( [3.3] ), except 
that the parameter domain now gets extended in the way indicated at the beginning 
of the previous section. 

In order to describe the evaluation formula for (pa>Pa)jv\ it is convenient to 
introduce dual parameters t r , r = 0, . . . ,3 that are related to the parameters t r , 
r = 0, . . . , 3 in the following way 

(tahtcUq- 1 ) 1 / 2 , 
(taht-H-'q) 1 / 2 , 
(tat^t^q) 1 / 2 , 
(tJ^t-Hdq) 1 / 2 , 



(4.1) 



ta 
tb 

i c 



where c and d denote the two indices that complement the indices a and b en- 
tering the truncation condition (3.6) such that {a,b,c,d} = {0,1,2,3} (cf. also 
the comment just after Lemma [3.1| ). It is worthwhile noticing that the parameter 
transformation defining the dual parameters in (4.1) is an involution (the duals of t r 
bring us back to t r , r = 0, . . . ,3) and, furthermore, that the dual parameters t r sat- 
isfy the truncation condition (3J3) if the parameters t r do so (because t a tb = t a tb)- 



We now form discrete Harish-Chandra-like c-functions C± R that are dual to the 



c-functions C± R ( 3.11 ), ( 3.12 ) (i.e., C±" is obtained from CJ H by replacing the 
parameters to, . . . by the dual parameters Iq, . . . , £3) 

i-l. 



(4.2) Cf(X) 



(4.3) C1 R {\) 



with 
(4.4) 



ao(A) J] I 

l<j<k<n 

>< n 

l<j<n 

c (A) n i 

l<j<fe<n 

>< n . 

l<j<n \ 

s (a)= n 

l<j<n 



(TjT k ;q)x 3 +\ k ( 



T j T k 



(ttjT k ;q)\ j+ \ k (tfjf k 1 \q)\ ] -\ k J 

Vllo<r<3(^3;^ 3 / 

(t- 1 qT j f k ; q)x J+ \ k (t^ 1 qT j T k r 1 ;q) Xj _ Xk 
(qfjfk;q)x j+ x k {qf j T^ 1 ;q) Xj -\ k 

v (qT 2 ;q) 2 \ 3 J 



<"^(wy 3 <r 1 ) 1/2 



and (cf. (|3J)) 
(4.5) 



t n -H a =t n -J(t t 1 t 2 t 3 q- 1 ) 1 / 2 . j = l, 



It is important to convince oneself that, despite the appearances of square roots in 
the definitions of l r (O) and tj (|J), the functions C^ H (A) (|J), (|J) (including 
the common factor cq (1-4)) are rational in the parameters t, to, . . . , <3 and q. The 
point is that in the above expressions the quantities t r and fj always occur in 



rational combinations of the form fjf k 



±i 



and that {tot^tsq' 



-l\l/2 



(iaibicidq 1 ) 1 ^ 2 = t a . In this sense the square roots in t r ( [4.l| ) have merely a 
formal meaning and these dual parameters were introduced mostly for notational 
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convenience and to emphasize the duality between the expressions for the above 
c-functions C*J R and the c-functions CJ R appearing in the previous section. 

Our main object of interest in this section will be a function A/" qR (A) on the cone 
A dJh defined by 



(4.6) 



7V qR (A) 



(7 qR (A) 
Cf(X)' 



It is clear from (the dual version of) Lemma |3.1| (part a.) and the above comments 
that C± R (A) and hence 7V qR (A) are well-defined nonzero rational expressions in 
the parameters t, to, . . . *3 and q for all A £ A. Furthermore, since i Q i(, = t a tb the 



truncation condition (3.6) reads the same in the dual parameters i r as in the original 
parameters t r . Thus, we may apply again (the dual version) of Lemma (part 
b.) to conclude that for A G An these expressions are well-defined and nonzero as 
rational expressions in the parameters t, to, . . A3 and q subject to the truncation 
condition (|3.6|). 

The following theorem provides a summation formula expressing {pxjPx}^ in 
terms of (1, l)^ (which corresponds to A = 0). 



Theorem 4.1. For parameters subject to the truncation condition (3.6) one has 
that 



(4.7) (p AjPA >Jf =JV« a (A)<l,l>J? 
withAf qR (\) given by ([D]) and 

(4.8) (1,1)??= E A? » 



(A e A N ), 



(where (■,■)% and lS. qR (y) are defined by (3.17) and (3.10), respectively) 



For an outline of the proof and further details regarding the Pieri type recurrence 
formulas lying at the basis of this proof the reader is referred to Appendix [b]. 

It is clear that the summation formula (4/7) should be thought of as a set of 
identities for the multivariab le A skey- Wilson polynomials with parameters subject 
to the t runcation condition (|3.6|) , which complements the orthogona lity identities 
111 lfT8|) . For n = 1 the formula reduces to the known (cf. |AWl| , |GRf l) summation 
formula 

(4.9) 

(A = 0,. 
(4.10) 



E 

0<v<N 



{px{t a q v )f/^{ V ) 



0?;?)aI1 



0<r<s<3 



{t r t s ;q) x 



(totit 2 t3q x 1 ;q)\ (£0*1*2*3; q)2\ 

N), with A qR (^) taken from ( |3.15| ) and 

E AqR (^)' 



(M>5? 



0<v<N 



for the monic g-Racah polynomial given by p\(z) (2.9) with parameters subject to 
the condition t a tb = q 



-N 



Although in principle the expression in (4.8) is in itself already explicit as a finite 



sum of explicitly given terms (this in contrast to the sum {pxyPx}^ with A ^ 0, 
for which in general one does not know the terms in a very explicit form unless 
n = 1), it would be interesting to also evaluate the sum (1, 1)° ^ in terms of a 
product formula. For n=l this was done by Askey and Wilson in [AW1] by means 
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of a g-Dougall summation formula for a very- well-poised (terminating) 6<fc series 
(see also |GR ), entailing 



(4.11) 



Q<v<N 



q;q) 



N 



(ibi c , ibid', q)N (t a tc 1 q, t a t d 1 q; q)N 



(where A qR (^) is given by ( p. 15 ) and t a tb = q~ N )- For arbitrary number of variables 
summing (1,1)^ amounts to the evaluation of a finite Selberg type q- Jackson inte- 
gral (related to the root system BC n ). Similar (but infinite) Selberg type q- Jackson 
integrals appear in the recent works of Aomoto and Ito jAo| , |] . More specifically, 
the q- Jackson integral corresponding to (1,1)^ constitutes a finite analogue of an 
infinite q- Jackson Selberg integral related to the root system BC n that belongs 
to the same class as those studied in [ |Ao| , ||. The connection with the q- Jackson 
integrals considered by Aomoto and Ito becomes particularly transparent after re- 
calling (see Remark 3.1) that for < \q\ < 1 the weights A qR (^) may be thought 
of as the restriction of a function A qR (z) of the form in ( |3.22 ) to the grid points 
rq", v G Ajv. 

We see from Theorem 4T that the quantities (j>\,p\)^ , A G An arc nonzero 
rational ex pre ssions in the parameters t, t$, . . . ,£3 and q (subject the truncation 
condition ( |3.6| )) because the factors A/" qR (A) are nonzero. (Clearly the common 
factor (1, l)f in (JO) is nonzero because otherwise the bilinear form (•, ■) N ^' would 
vanish identically on the space Jiff .) This checks with Proposition 3.4, which stated 
that for generic parameters (subject to the truncation condition) the bilinear form 
{•,•)% is nondegenerate on the space fC^. 

The orthonormalization constants are given by square roots of (p\,p\)'^ i , A G 
Aj\r. After division by these constants the basis p\, A G An becomes an orthonormal 
basis for the space Tiff with respect to the bilinear form (•,•)??• The fact that 
(generically) the orthogonal basis pj, A 6 A^ can be turned into an orthonormal 
basis for Tiff (or, equivalently, that the quantities (p\,p\)]^' , A G A^r are nonzero) 
implies that (generically) no linear dependences arise between the functions p\, \ <E 
An when being restricted to the grid rq v , v G An- This is of course precisely what 
we already saw in Lemma 3.5, and what — because of dimensional considerations — 
implied that (generically) any function over the grid could be represented exactly 
by a unique function in Tiff (Proposition 3.6). 

We also saw already in Sect ion |3| that for parameters subject to the truncation 
condition the weights A qR (^) (3.10) vanish when v € A \ A^v, because then the o 
function C^iy) (3.11) became infinite. Since the truncation condition ( |3.6| ) reads 
the same in the dual parameters t r (4.1) as in the original parameters t r (recal 



igib — t a tb), the c-function C+ R (A) (|4.2|) also becomes infinite for parameters subject 



to the truncation condition when A G A \ An- Hence, the sum (p\,p\)n^ (4.7) 
vanishes if the parameters satisfy the truncation condition (3.6) when A G A \ A^v- 
The vanishing of the sum (pA,PA)5f in this situation should of course not come as 
a surprise in view of the previously noted fact that the polynomials in question are 
actually zero on the grid points Tq v , v G An for these parameters (Proposition 3.7). 

It is possible to capture the orthogonality relations (3.18) together with the 
orthonormalization formulas (4.7) in terms of a purely linear-algebraic formulation. 
In this language Theorem 3.2 and Theorem 4.1 boil down to the property of the 
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matrix 

/- \i/ 2 / \V 2 
(A^^)) ( A qR (^) ] 

(4.12) [K^ veAN with ^sP^P 7 



;i, i) qR 



1/2 



being an orthogonal matrix for generic parameters subject to the truncation con- 
dition (3.6). Here we have used the renormalized polynomial P^{z) = C"^ R (fi)p u (z) 



(cf. Appendix |b|) and the weights for the 'Plancherel' measure (cf. ( |3.10[) ) 

A qR (u) = - 

and, we have again (cf. ( [3. 20] )) ordered the rows and columns of the matrix accord- 



ing to an arbitrary total extension of the partial order in (2.5 ). T he orthogonality 



(and thus invertibility) of the matrix K = [K^^^^^Ax in d4.12|) checks with the 



previously noted invertibility of the matrix M given by (3.20) because — as should 
be clear from our discussion — one has that K = LMD with L triangular and D 
diagonal and both L and D (generically) nonsingular (as the matrix elements on 
their respective diagonals are nonzero as meromorphic functions in the parameters 
subject to the truncation condition). 

Remark 4.1. The normalization of the polynomials P^{z) — C'f < '(ii)p fl (z) is such 
that P^(t) — 1 (see Appendix [b]). Moreover, the renormalized multivariable Askey- 
Wilson polynomials P^(z) satisfy Macdonald's duality property (see also Appen- 
dix |b|) stating that 

(4.13) P^rq") = P„(f^), 

where P u {z) = C+ R (^)p„(z) denotes the renormalized multivariable Askey- Wilson 
polynomial dual to P v {z) = C^ R '(i')p u (z), i.e. with the parameters t r being replaced 
by the dual parameters i r ( |4.l| ). (For v — the equality in ( 4.13| ) reduces to the 



evaluation formula P^ij) = 1.) The duality relation ( 4.13j ) for the polynomials 



gives rise to a similar duality property for the matrix [-K^J^eAjv (4-12), viz. 
transposition of [K^^]^^(za n leads one to the matrix [-K^J^eAjv m which the 
parameters t r are replaced by the dual parameters t r . Here we also used that 
the orthogonality of the matrices fixes their normalization uniquely. In particular 
one sees that as a consequence the sum (1,1)^ is invariant with respect to the 



transformation t r — > t r (for parameters satisfying the truncation condition (3.6)). 
For n = 1 this property of (1,1)^ boils down to the equality 

U 14 ) (i 2 a q,i-% 1 q;q) N = (tlq,t- 1 t^ 1 q;q) N 

(i a ic 1 q,i a i d 1 q\q)N {tJc l q, t a t^ x q; <?)jv ' 

which is not difficult to check directly using ( |3.6| ), fl4.l| ) and the transformation 
property for q-shifted factorials 

(a 1 q~ N ;q) N /(a 2 q~ N ;q)N = (ai/a 2 ) N (a~ 1 q; q) N /(a~ 1 q; q) N . 

In the special case that the parameters satisfy the additional constraint t a q = 
tbt c td one has that t r = t r (r = 0,... ,3). Hence, we are then in a self-dual 



situation Pp,(z) — P(t(z) and the matrix [K^ V \^ V £A N (4-12) — in addition to being 
orthogonal — now also becomes symmetric. 
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5. Transition to Racah type polynomials 



We will now study the transition from the basic hypergeometric level to the 
hypergeometric level. To this end we substitute the variables 

(5.1) Zj = q x >, j = l,...,n 

and in addition perform a reparametrization of the form 



(5.2) t = q g , t r = q a % r = 0,...,3. 

After these substitutions and division by a constant factor (1 — q) 2 , the g-difference 
operator D (^J) passes for q — > 1 over into a second order difference operator given 

by 



(5.3) 
where 



l<j<n 



(2exj) (1 + 2exj) 

>< n 

l<fc<r, 



Xk 



Xk 



g + ex,j 



Xk 



Xk 



±1 



and with the action of the operators Tj, j = 1, . . . , n being of the form 

\-^jf) \X\t • • • ) X n ) — f\Xl ) • • • ) Xj— x , Xj -\- 1, Xj-j-i, . . . , x n ) . 



It turns out (sec |D3|) that the difference operator D (5.3) is triangular with 
respect to the (partially ordered) basis of symmetrized monomials for the space 



C 4 



'1> • 



(consisting of the permutation-invariant and even polynomials in 



the variables x\, 



(5.4) 

where 
(5.5) 



Drh\ = 



,x n ). Specifically, one has that 

Ex^rh^ with E\^ G C[g, go, gi, ff2, 53] 



E 



fh\{x) = E ■'' 
^es„(A) 



2/<i 

1 



r 2^ 



A E A. 



Here the summation in (5__5) is meant over the orbit of A 6 A (2J2) under the action 
of the permutation group S n (which permutes the vector components Ai, . . . , X n ) 
and the partial order of the cone A (2.2) is taken to be the same as before (see (2.3)). 
The d iago nal matrix elements E \.\ in (5.4) (which can be obtained for q — * 1 from 
E\^\ (2J:) after substitution of ( |5.2| ) and division by (1 — q) 2 ) read explicitly 



(5.6) 



l<j<n 



A e A 



with 

(5-7) = (n-j)g + (g Q + g 1 +g 2 +g3-l)/2, j = l,...,n. 

The triangularity of the difference operator again reduces the corresponding eigen- 
value problem in the space of the permutation-invariant and even polynomials to an 
in essence finite-dimensional problem. Although the eigenvalues (^1]) are no longer 
nondegenerate, it still remains true that E\ \ ^ E^^ as polynomial expression in 
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the parameters g, go, . . . , g$, if A ^ fi and A, fx are comparable with respect to the 
partial order ( |2.3| ). Fortunately, this is already sufficient to single out the eigenfunc- 
tions uniquely by means of conditions analogous to those entering the definition of 
the multivariable Askey- Wilson polynomials p\ in Section 0. It turns out (see [ D3 1 ) 
that in the present case we are in fact dealing with a multivariable analogue of the 
Wilson polynomials |W|. 



Definition. The multivariable Wilson polynomial associated with a (dominant 
weight) vector A 6 A (2.2) is the (unique) monic permutation invariant and even 
polynomial of the form 



(5.8a) px(x) = fh\(x) + 

such that 
(5.8b) 



Dpx = Ex,xpx- 



We can again represent these polynomials in terms of a formula of the type (2 



(5.9) 



Px 



n 



D — E, 



Ex, 



E„ 



mx- 



Furthermore, one may always replace the monomial basis in the r.h.s. of such a 
formula by any other basis related to it via a unitriangular transformation (this is 
immediate from the comments following ( |2.8|) that proved the validity of this type 
of representations for the polynomials). In [D3 it was demonstrated, by performing 
a suitable unitriangular transformation of the basis elements, that formula (2.8) for 



the multivariable Askey- Wilson polynomials tends to formula (5.9) for the multi- 
variable Wilson polynomials in the limit q — > 1 after substitution of ( |5.l| ), (5.2) and 
division by a constant factor (1 — q) 2 ' A L Here we have used the (standard) notation 

(5.10) |A| = Ai + --- + A„. 



(The point of the unitriangular transformation is that the basis elements (2.1) all 
collapse into constant functions for q — ► 1 after the substitution (5.1); by temporar- 
ily passing to a basis of TL W with elements of the form 



(5.11) 



E 

/i£Sr l (A) 



(zi + z^ 1 - 2)^ • • • (z n + z- 1 - 2Y 



A g A, 



it is seen that a nontrivial limit is obtained after one divides out the constant factor 
(1-q) 2 ^.) 

The upshot is that there exists the following limiting relation between the mul- 
tivariable Askey- Wilson polynomials px of Section and their Wilson type coun- 



terparts px of the present section [D3|. 

Proposition 5.1. For Askey-Wilson parameters given by ( |5.2| ) ; one has that 

(5.12) px(x) = ]im(l-q)- 2 ^px(q x ) AgA 

g-fl 

(where q x = (q Xl , . . . ,q Xn ))- 

In |D3| the orthogonality properties of polynomials px (x) were investigated with 
respect to a continuous Wilson type weight function A w . For parameters satisfying 

(5.13) g>0, Bje(g r ) > (r = 0,1,2,3), 



IS 



J. F. VAN DIEJEN AND J. V. STOKMAN 



with possible non-real parameters g r occurring in complex conjugate pairs, the 
relevant orthogonality relations read 

/oo />oo 
■■• / p\(ix)fi[_ l (ix) A w (x) dxi ■ ■ ■ dx n = if 
-00 J — 00 

where 

A w ( ) = TT r (9 + i{£iXj +£2Xk)) 
{ > J} k<n r(i(s 1 x j +e 2 x k )) 

El ,E2=ll 



n 



T(g + isxj)T(gi + iexj)r(g 2 + iexj)T(g 3 + iexj) 



Y(2iex i 

l<j<n V J 

e=±l 



(with r(-) denoting the gamma function). 

We will now apply the limiting relation of Proposition 5.1 to the results of Sec- 
tion^ and ^ to infer that for generic parameters subject to the truncations condition 

(5.15) {n-l)g + g a + g b + N = 

(with N a nonnegati ve in teger and a, b £ {0, 1, 2, 3} such that a ^ b) the polyno- 
mials p\{x), A 6 An (3J5) constitute an orthogonal basis for the finite-dimensional 
space 

(5.16) H% = Span{mA}AeA N 
endowed with a nondegenerate bilinear form determined by 



(5.17) (f,g)% = J2 f(p + v)g(p + v)& R (v) (f,gen 



Here the vector p is of the form (c.f. also the vector p in Remark 3.2 with r given 
by ( |3.8| ) and parameters taken from (5.2)) 

(5.18) p = (px, . . . ,p n ) with Pj = (n - j)g + g a , j = l,...,n 
and the weights are given by 

(5 - 19) A » = cmcwy 

with 

(5.20) Cliv) = TT ( (Pj+Pk)u ] +, k (Po-Pk) Vj -v k \ 

!</<!<„ + Pi + W + Pi - Pk)vj-Vk J 

&Pj)2v 3 \ 



>< n 

l<j<n 



(5.21) C R (v) = n 



rio<r<3(fr + Pj) Vj ) ' 
0--9 + Pj + Pk) Vj + Vk (1 - g + Pj - Pk) Vj -v 



l<j <fe <„ x (! + P3 + Pk)u ]+ y k (1 + Pj ~ Pk) V] -v k 



>< n 

l<j<n 



no<r<3( 1 -.9r+Pj)^ 



(1 + 2^)2^. 

where we have used Pochhammer symbols defined by 
(eti, . . . , a s ) m = (ai) m • • • (a s ) m , 

(a) m = a(a + 1) • • • (a + m — 1) 
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(with (a)o = 1)- 

In order to describe the corresponding orthonormalization constants we again 
need dual parameters 

' 1 1 1 1 
1 1-1-1 
1-1 1-1 
1-1-1 1 

({a, b, c, d} = {0, 1, 2, 3}) and a function 7V R (A) on the cone A (|2.2|) 

C* R (A) 



(5.22) 



(9a\ 


( 


9b 


1 


9c 


~ 2 


W 


\ 



\ 


(9a\ 








9b 


1 


1 




9c 


+ 2 


1 


I 


\W 




V 1 / 



(5.23) 



Af R (\) = 



C?(A) 



that is governed by c-functions dual to C R (5.20), ( ]5.21 ) 
(5 24)(7 R (A)- 77 ^ (A? + P&K+A* (pj ~ Pk)\ 3 -\ k 




(5.25) C* R (A) 



n 

l<j<k<n 



(.9 + A, + /5fe)A J+ A fc (g + Pj - Pk)\ 3 -\ k 

]lo<r<3(.9r +Pj)a 3 y ' 
(1 - .9 + Pj + Pfe)A 3 +A fc (1 - .9 + Pj - Pfc)A 3 -A fc 



>< n 

l<j<« 



(1 + Pj +/5fe)A 3 +A fc (1 + Pj 



Pk)\ 3 -\ k 



(1 + 2ft) 2 A 3 



(Recall that the components of the vector p are given by (5.7), so pj — (n—j)g+g a ). 

The following theorem, which describes the orthogonality properties of the poly- 
nomials px, A G An with respect to the bilinear form (■, ■ ) R / , is an immediate 



consequence of the application of Proposition 5.1 to Theorem 3.2 and Theorem 4.1 



Theorem 5.2. For parameters subject to the truncation condition (5.15) one has 
that 



(5.26) 
and that 
(5.27) 



{Px,Px) R N =M R (\) (1,1)% 



(X,p 6 A N ) 
(A e A w ) 



with M H (\) given by (5.23) and 

(5.28) <M>& = 



E 

vGAjv 



(where (■,■)% and A R (v) are defined by (5.17) and (5.1£), respectively). 

(The formulas (|J|) and (|5.27|) should again be interpreted as equalities between 
expressions that are rational in parameters subject to the truncation condition.) 
To verify the theorem it suffices to infer that for q — > 1 and parameters given 

(multiplied by a factor 



by Q the c-functions Cf-(v) ( pUlp and C^ R (A) 
(1 - q) 2 ^ and (1 - q) 2 ^ , respectively) converge to C^(v) (|5~2p|) and C*(A) (^24|) ; 
and that, similarly, the c-functions C"i R (z/) ( |3.12| ) and C"i R (A) (Q (divided by a 
factor (l-q) 2|l/| and (1 -g) 2 ^ , respectively) tend to C R (v) ( [5~2l"l ) and C* R (A) fl5l>5l) 
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in this limit. To this end one simply uses that the (renormalized) g-shifted factorial 
(a; q) m / (1 — q) m converges to the Pochhammer symbol (a) m when q tends to one. 
Furthermore, for the parameters ( |5.2| ) the truncation condition (3.6) amounts to 
the condition (n — l)g + g a + gb + N = (mod 2ni/ log(g)), which entails ( 5.15 ) in 
the limit q — > 1. 



It is important to again convince oneself that the formulas (5.26) and (5.27) 
are indeed well-defined as rational expressions in the parameters g and go, ■ ■ ■ , <?3 
subject to the truncation condition ( 5.1 5| ) (i.e., no denominator becomes identical to 
zero) and, furthermore, that the r.h.s. of ( |5.27 ) is nonzero as a rational expression 



in these parameters. We thus have that the bilinear form (•, -) N is nondegenerate on 



the space TL N for generic parameters subject to the truncation condition (5.15) and 



that any function defined on the grid points p+A, A £ A at can be represented exactly 
by the restriction to the grid of a unique function in the space H N . (The restriction 
to the grid points of the orthonormalized basis ((pa,Pa) R /)~ 1//2 Pa(2;), A £ An for 
Ti 1 ^ yields an orthonormal basis for the space of functions over the grid p + A at 
endowed with the nondegenerate bilinear form X^eAw f(P V )9(P i / )A r (j/).) 
Finally, as degeneration of Proposition 3.7 we arrive at a similar statement for the 
polynomials p\ . 



Proposition 5.3. For parameters subject to the truncation condition ( 5.15 ) one 
has that 



(5.29) 

ifX £ A\ A N . 



p\ (p + v) =0 for v £ A 



A 



The corresponding linear-algebraic formulation of Theorem 5.2 states that the 
matrix (cf. fl4.12j )) 

/' \V2/ \l/2 

(a r (m)) (a») 

[VUa, with = P^Tq")^ 



(5.30) 



((1,1) 



1/2 



is orthogonal for generic parameters g, go, ... ,173 subject to the truncation con- 
dition ( 5.15 ). Here P M is the renormalized polynomial Pp,(x) = C^:(p) p fJl (x) and 
A R denotes the weight function for the discrete 'Plancherel' measure A R (/i) = 



l/(Cf(p) C R (p)) dual to the weight function A R (|5T9| ). 

As a side remark we mention that the orthogonality of [K^^n^^Kn implies that 
the matrix 



(5.31) 



M = [m M (p + ^)] Mjl / G A r 



has a determinant that does not vanish as a polynomial in the parameters g and 
9a- This is because M is related to K = [-^^,i/]^,i/eAjv ( M0|) by K = LMD, 
where L is triangular and D is diagonal with both matrices having elements on 
the diagonal that do not vanish as meromorphic functions in the parameters g and 
go, ■ ■ ■ ,33. It is interesting to observe that the invertibility of the matrix M does 
not seem so easily established directly (i.e. without using K), as was the case 
when dealing with its q- version M — [m ll {Tq v )] tJli ^ e x N ( 3.20| ) just after Lemma 3^ 



(The problem is of course that here in the degenerate version we have lost the 
possibility to play with the parameter q.) Notice, however, that in the special case 



of only one single variable the matrix M (5.31) becomes a Vandermonde matrix 
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[(g a + v) 2fJ "]o<fj,.v<N , from which it immediately follows that the determinant is 
nonzero for generic parameter values. 

The difference equation ( 5.8b| ) tells us that for n — 1 the polynomials p\ reduce 
to the monic Wilson polynomials |w| 

,„„v - / \ (90 +91, 9o +92,90 +gs)x „ 

(5.32) px(x) = 



(90 
4F3 



9i + 92 + 93 + A - 1) A 
-A, .go + 9i + 92 + 93 + A - 1, g a + 
9o + 9i, 9o + 92, 9o + 93 



9o 



where we have used standard notation for the hypergeometric series (see e.g. [ AW2 

grJ) 



TP 



CLl, ■ 

bx,. 



,a r 
.6, 



;z = 



(ai, . . . ,a r ) k z k 
,b s )k k\ 



^ (61 

k=0 v 



(The explicit hypergeometric representation ( |5.32 ) for p~\ in the case of one variable 
also follows from Proposition 5.1 and the corresponding basic hypergeometric for- 
mula for px in (2.9).) The identities in Theorem 5.2 amount in this special case to 
the discrete orthogonality relations for the monic Wilson polynomials ( |5.32 ) subject 
to the parameter condition g a + gi + N = [|w| 

is)A B (v) =0 for X^fi 



E 

0<v<N 



p\(g a + v)Pp,{g a 



(A,/xe{0,. 



E 

0<u<N 



, N}) and 
P\(9a + v)p\{g a 



-v)A*-(v) = 

Xl Uo<r<s<3(9r+ 9s)x 



(A G {0, . . 

<M> R 



JY 



(go + 9i +92 + 93 + A 
, AT}), where 

= E A »> 



1)A (50 + gi + .92 + ff3)2A 



0<v<N 

- d + - 

g a 



(go + ga,gi + ga, 32 + ga, <?3 + ga)i 



(I - go + g a , 1 - gi + g a , 1 - .92 + .g Q , l - 5.3 + flh)* 
In the finite-dimensional case with the truncation condition g a + gj, + N = and 



discrete orthogonality properties, the one- variable polynomials f>\ in (5.32) are usu- 
ally referred to as Racah polynomials rather than Wilson polynomials. The nor- 
malization factor (1, for the one- variable Racah polynomials can be evaluated 
in product form by means of a summation formula for a very-well-poiscd (termi- 
nating) 5F4 series due to Dougall, which entails |w| (see also GR ] for the Dougall 
5-F4 summation formula) 

(l + 2g a ,l-g c -gd)N 



(5.33) 



E A » 



(1 



9c 1 + g a - go) 



N 



Remark 5.1. The duality r elat ions for the (renormalized) multivariable Askey- Wils- 
on polynomials in Remark |4 - 1| give in the limit q — > 1 rise to analogous duality rela- 
tions for the (renormalized) multivariable Wilson polynomials P^(x) = C+(f/,)pn(x) 
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(5.34) P ll (p + v) = P v tf + n), 



where P v {x) — C^{v)p v {x) is the dual of P v { x) — C^:(v)p v (x) with the parameters 
g r being replaced by the dual parameters g r ( 5.22| ). (For v = the duality relation 



( 5.34) reduces to the evaluation formula P^ip) = 1 characterizing the normalization 
of the polynomials P fJ- (z).) Just as in the q-case, the duality properties for the 



polynomials are again inherited by the matrix \K u , u ] u , ue / i _ N in (5.30). We now have 
(for parameters satisfying the truncation condition (|5.15|) ) that transposition of 
the matrix [-F^i/J^gAn amounts to the parameter transformation g r — > g r and, in 
particular, that the sum (1, 1)^ is invariant with respect to such a transformation 
of the parameters. For parameters satisfying the additional constraint g a — gt — 
9c — 9d = — 1 we have that g r = g r . Hence, in that case we are again in a self- 
dual situation (i.e. P v (x) = P v {x)) and the matrix [K fi ^] fJ _^ e \ N now also becomes 
symmetric (in addition to being orthogonal). 



Remark 5.2. It turns out that our multivariable Racah polynomials (i.e., the mul- 
t ivari able polynomials p\(x) with parameters subject to the truncation condition 
( 5.15| )) are not the first generalization of Wilson's one- variable Racah polynomials 



to the case of several variables. Already several years ago Gustafson reported on 
a finite system of multivariable orthogonal polynomials with discrete orthogonal- 
ity measure tied to the so-called multiplicity-free Racah coefficients for the group 



U(m + 1) [Gu]. In the rank one situation (m = 1), the orthogonal polynomials in 
question can be reduced to the Racah polynomials of (jw|. It is not clear (at least 
not to us), though an interesting question, whether also for higher rank Gustafson's 
multivariable Racah polynomials may be linked with the multivariable Racah poly- 
nomials of the present paper. From the explicit expressions for the weight functions 
it seems that both approaches generalize the one- variable Racah polynomials to sev- 
eral variables along very different directions. However, at present we are not able 
to rule out completely the possibility that there might not be some transformation 
connecting the two approaches. Should there indeed exist such a connection (which, 
however, would seem more likely after trading the group U(m + 1) for Sp(m) say), 
then this would imply an interesting link between Gustafson's group-theoretical 
program and a degenerate case of the Macdonald theory. 



Remark 5.3. The reader may find it illuminating to view the multivariable q-Racah 
polynomials p\ in Sections [3] and | as a trigonometric version of the multivariable 
Racah polynomials p\ in the present section. If we substitute Askey- Wilson pa- 
rameters in accordance with (|5.2|) and set q — e la , then we may rewrite the weight 



function A qR ( |3.10D as 

A qR (i 



Cf(y) C^{vY 
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with 

CT» = (-4)- 




and 

C q V) = (-4)1"! 



x 



n 



l<j<k<n 

>< n 

l<j<n 



(1 - g + pj + p k : gin. a ) Vj+ „ k (1 - g + pj - p k : sin a )^._ l 
(1 + pj + pk : sina)^^ (1 + pj - p k ■■ sin a )^_^ 

n <r<3( 1 - 9r + Pj ■ Sm Q )„. 



(1 + 2pj : sia a ) 2vj 

(where the components of the vector p are given by ( 5. IS )). In the above formula 
we have used 'trigonometric Pochhammer symbols' defined by 

(5.35) (a : sin a ) m = sin Q (a) sin a (a + 1) • • • sin Q (a + m — 1) 

with (a : sin Q )o = 1 and sin a (£) = sin(a£/2). Furthermore, we also arrive 
at corresponding trigonometric expressions for A qR (A) = l/(C , ^ R (A)C'!i R (A)) and 
7V qR (A) = C^\X)/Cf{\), which are governed by dual c-functions^jf obtained 



by replacing the parameters g r by g r ( 5.22 ) and the vector p by p (p.7|) 



It is manifest from these representations that A qR , A qR and jV qn can be in- 
terpreted as trigonometric versions of A R (|5l9|) , A R (see ( pc|) ) and 7V R ( ^23|) , 



respectively. The transition q — > 1 corresponds to the limit a — > in which the 
renormalized trigonometric Pochhammer symbols (2/a) m (a : sin a ) m go over in the 
ordinary Pochhammer symbols (a) m . As far as the polynomials are concerned, we 
see that p\(e lax ) becomes a trigonometric polynomial in the variables x%, ... , x n 
with period 2ir/a. Proposition describes the rational limit in which the period 



of the trigonometric functions tends to infinity (cf. [D3]) 



Notice also that the modified monomial basis elements in (5.11) are in the trigono- 
metric notation of the form (— 4)' A 'm,\(sin Q (2;)) and converge, after division by the 
constants («a) 2 ' A ', to fh\(x) for a — » . Finally, in the trigonometric coordinates 
the grid points on which the discrete orthogonality measure A qR for the polyno- 
mials p\{e lax ) is supported become of the form p + u, v € Ajy- This means, in 
particular, that in these coordinates the grid points do not move when performing 
the limit q — > 1 (or cquivalcntly a — > 0). 

6. POSITIVITY DOMAIN FOR \q\ = 1 AND THE q-RACAH TRANSFORM 

In the preceding sections we have viewed the polynomials and other objects of 
interest (such as the discrete weight function and the difference equation) as rational 
expressions in the parameters. As a consequence, we arrived at their properties for 
generic values of the parameters. It is clear, however, that the above generic picture 
does not hold for all values of the parameters. A most drastic way in which the 
generic picture breaks down occurs when the cardinality of the grid rq", v 6 A^ 
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becomes less than the dimension of the space Tiff (3.16) (or, equivalently, when 



it becomes less than the number of points in the alcove Ajv ( p.5[ )). This may for 
instance happen when q M = 1 for M E {1, ... , N}. The special case M = 1 (so 
q = 1) corresponds of course precisely to what we already analyzed in more detail in 
the previous section by means of a limit transition. In this section, however, rather 
than to present any in-depth analysis of the case that q is such a root of unity 



(cf. [C2, Ki, SpZ]), we will focus on a parameter domain with \q\ = 1 for which 
the generic picture sketched in the preceding sections does apply for all parameter 
values in the domain and, moreover, for which the weig ht function A qR (i/) fl3.10| ) 
becomes real- valued and positive when v lies in the alcove Ajy- 

To the describe the positivity domain for \q\ — 1 it is convenient to perform the 
following trigonometric substitution of the parameters (cf. Remark |5.3|) 



(6.1) 



q = e' u 
to = e 1 



tb 



= e lag , 

— _p»«g6 



tc _ e ia( Sc +l/2) td _ _ e ia(gd + l/2) _ 



The weights A qR ( 3.10| ) can then be rewritten with the aid of the trigonometric 
Pochhammer symbols (cf. Remark 5.3) 



(6.2a) (<zi, . . .a s : sin Q ) m = (cti : sin Q ) m • • • (a s : sin Q ) m 

(6.2b) (a : sin Q ) m = sin Q (a) sin Q (a + 1) • ■ ■ sin Q (a + m — 1) 

(6.2c) (ai, . . . a s : cos Q ) m = (a x : cos a ) m • • • (a„ : cos a ) m 

(6. 2d) (a : cos Q ) m = coSo,(a) cos a (a + 1) • • • cos Q (a + m — 1) 

where (a : sin Q ) = (a : cos a ) = 1 and sin a (£) = sin(a^/2), cos Q (^) = cos(a^/2). 
Specifically, we have 



(6.3) 
with 

c?(»)= n 



A qR (^) 



1 



><n 

l<j<n 



l<j<k<n 



(Pj + Pk ■ sm a )„ j+ „ k (pj - p k : sinq)^.-^ 
(.9 + Pj + Pk ■ sin a (g + pj - p k : sm a ) Vj -u k 

(pj, 1/2 + pj : siriq)^ (pj, 1/2 + pj : cos a )^. 

(ga + Pj,Sc + 1/2 + Pj ■■ sin a )^(g b + pj,g d + 1/2 + pj : co 

(1 - g + pj + p k : sin a )^ +!/Ji (l-g+ pj - p k : sitia.)^-^ 



l<j<k<n 



><n 

l<j<n 

and 



(1 + pj + p k : sino,)^^ (1 + pj - pk : sin Q )^._, yfc 

(1 - g a + Pj, 1/2 - g c + Pj ■ shiq)^ (1 - g fc + pj, 1/2 - g d + pj : cos a 
(1 + pj, 1/2 + pj : sm a )^(l + pj, 1/2 + pj : cos a ) Vj 

Pj = (n-j)g + g a - 



The corresponding dual objects CI 11 and are again obtained by replacing the 
parameters g r by dual parameters g r with 

(6.4) 







f 1 


1 


1 


1 






fga\ 


Sb 


1 


1 


1 


-1 


-1 






gb 


gc 


~ 2 


1 


-1 


1 


-1 






gc 


\gdj 




V 1 


-1 


-1 


1 


J 




{gdj 
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and the vector p by p with pj = (n — j)g + g a 



Proposition 6.1. For parameters satisfying the constraints 

7T 

(6.5a) a > 0, .g > 0, 0<g a ,g 6 <-, -g a < g c < g a 

a 

and £/ie truncation condition 



-gb <gd< gb 



(6.5b) 

one has that 



(n-l)g + g a + g b + N = -, 

a 

< Cf{u) < oo for v G A^v, 
< CjfO) < oo /or /i G A N . 

Proof. Notice that for N = the proposition is valid trivially, because in that case 
CJ R ,C± R = 1 (and An = {0}). Let us from now on assume that N is positive 
and, furthermore, let us also temporarily assume that the parame ters q and g a , g& 
are nonzero. It is not very difficult to verify that the conditions ( |6.5a ) and ( 6.5b ) 
then imply that the arguments of the sinus functions in C± R (^) (with v G A^v) 
lie in the open interval (0,7r) and, similarly, that the arguments of the cosinus 
functions lie in the open interval (— 7r/2, 7r/2). Hence, it follows that the c-functions 
CJ R (^), v G Ajv are positive and finite for these parameters. When one or more 
of the parameters g,g a ,gb (or g c ,gd for that matter) become zero, one may cancel 
the factors in the numerator/denominator carrying that parameter (with the value 
zero) against the corresponding term in the denominator/numerator (thus resulting 
in a trivial unit factor). By proceeding in this manner one readily infers that the 
positivity and finiteness of the c-functions is conserved also when one or more of 
the parameters g and g a ,gb are allowed to become zero. 

The dual statement that C^'(fx) is positive for p, G Ajv is now immediate from 
the observation that the conditions ( |6.5a ), ( ]6.5b| ) are self-dual in the sense that they 
read the same in the parameters g a , gb, g c , gd as in the dual parameters g a , gb, g c , gd- 
Specifically, the conditions Q6.5a ), (6.5b) (with N > 0) are equivalent to the condi- 
tions 



(6.6a) 

and 

(6.6b) 



a > 0, g > 0, 



< ga,g6 < -, 

a 



< 



< 



-gb < gd < gb 



(n-l) 5 + g a + g fc + iV= -. 

a 



That the conditions ( 3.5a ) and ( 6.5b ) imply the dual condition s (3.6a ) and ( 6.6b ) 
is seen with the aid of the definition of the dual parameters in (6.4); that both the 
conditions and the dual conditions are actually equivalent is then clear from the 
fact that the parameter transformation determining the dual parameters in (|6.4[ ) is 
an involution. □ 



Corollary 6.2. For parameters subject to the conditions in Proposition 6.1, the 
weight functions A qR (u) and A qR (p) are positive and finite when v,p lie in the 
alcove An (3.5). 



The positivity of Proposition 6.1 implies that p\, A G An is well-defined not just 
generically but for all parameter values in the domain determined by the conditions 
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Proposition 6.3. The polynomials p\, A £ An are well-defined for all parameters 
( f3~l| ) with values in the domain det ermi ned by the conditions ( 5.5a ) and ( |6.5b| ) (i.e., 
the expansion coefficients ca, m in (2.7a) are regular for these parameter values). 



Proof. Let us first assume that the parameters are generic (complex say) but subject 
to the truncation condition ( |3.6[ ). For arbitrary A £ A^r, the bilinear form (■, 
restricted to the subspace Span{m M | \i < A} is nondegenerate as a consequence 
of Theorem 3.2 and Theorem 



4.1 



Using the orthogonality of the polynomials 
and the nondegeneracy of the bilinear form it is seen that one can characterize 
the multivariable g-Racah polynomial corresponding to a dominant weight vector 
A £ An as the unique polynomial of the form 

P\{z) = m x (z) + ^2 

/jGAjv, £i<A 

such that 



.(*) 



(P\,m 



/j./'n — 



for ji < A. 



In other words, p\ consists of m\ minus its (unique) orthogonal projection with 
respect to the (nondegenerate) bilinear form onto Span{m M | // < A}. We 

thus have the following inductive Gram-Schmidt-like formula for the (orthogonal) 
polynomials p\, A £ 



(6.7) 



px(z) = m\(z) - ^2 



(At this point it is helpful again to view Equation (6.7) as a rational identity in the 
parameters q, t and t r subject to the truncation condition (3.6).) Using induction 
on the weight A one sees from this inductive formula and Theorem h.l\ that the 
polynomial p\ is regular at parameter values where the c- functions , C+* are 
regular and nonzero. The proposition is then immediate from Proposition |6.l| . □ 



With the aid of Proposition 6T and 3J? it is seen that the results of Section Hand 
|] hold for all parameter values in the domain determined by the conditions (3.5a) 
and ( 3.5b| ). At some points the positivity of the measure will enable us to even 
formulate a somewhat stronger version of the results stated there. For instance, 
for the parameters in the positivity domain ( |6.5a ), (6.5b) the polynomial p\(e tax ) 
(x £ R") and thus also the matrix [K^^]^, v£ ^ N ( 4.12 ) is real. Hence, in addition 
to being orthogonal the matrix [K^^]^^^^ now also becomes unitary. (Th e rea l- 
valuedness of polynomials p\(e tax ) for parameters in the positivity domain ( |6.5a ), 
(3.5b) follows e.g. from the inductive formula (6.7) together with the observation 
that the (even) monomials m\(e lax ) are real.) The fact that both polynomials 
and weight function are real for parameters in the positivity domain ( |6.5a ), (6.5b) 
allows us to restrict the bilinear form (■, Ojv" Q3-17 ) to a real form (by restricting to 
the real vector space spanned by m\, A £ A at), which in turn can be extended to a 
positive definite sesquilinear form on the complex vector space ( ft.!6| ) (in the 
standard way). 

In the remainder of this section we will interpret the orthogonality and orthonor- 
malization properties of the multivariable q-Racah polynomials with parameters in 
the positivity domain in terms of a finite-dimensional discrete integral transforma- 
tion for grid functions. To this end we need to introduce some further notation. Let 
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L 2 (p + An, A qR ) be the finite-dimensional Hilbert space of complex functions over 
the grid points p + v, v G An endowed with the standard inner product determined 
by the positive weights A qR (i/), v G An 



(6-8) (/, g) A = V f(p + v) g(p + „) A qR (j/). 



E 



Similarly, the space L 2 (p + An, A qR ) denotes the corresponding dual Hilbert space 
consisting of the complex functions over the grid points p + p, p G An endowed 
with the standard inner product (•, ■) ^ determined by the positive weights A qR (/z), 
p G A N . We define the operator AC : L 2 (p + A N , A qR ) — ► L 2 (p + A N , A qR ) as the 
map with kernel 



(6.9) JC(p + n, p + v) 



_ P M (e M ( p+ ^)A qR (i/) 



i.e., the operator AC acts on a grid function / : p + An — * C as 
(6.10) (Kf)(p + fx)= lC{p + l x,p + v)f{p + v) 



(thus producing a function (AC/) : p + An — > C). Here we have employed the renor- 
malized multivariable Askey-Wilson/g-Racah polynomials P\(z) = C^ i (X)px{z) 



satisfying the normalization condition P\(t) = 1 (cf. Remark 4.1). We also define 
the dual map AC : L 2 (p + A N , A qR ) — ► L 2 (p + A N , A qR ) determined by the kernel 

(6.1D t iP+ ^ + ,^i^^m^ 

(P\{z) = (\)p\(z)) and acting on a function / : p + An — » C over the dual grid 
as 

(6.12) [tf){p + v)= AC(p + i/,p + /x)/(p + M ). 

The following theorem describes a discrete integral transform — the 'q-Racah 
transform' — between grid functions in L 2 (p + An, A qR ) and L 2 (p + An, A qR ) to- 
gether with its inversion formula. 



Theorem 6.4. For parameters subject to the conditions ( 5.5a ), ( 6.5b ) the map 
AC : L 2 {p + An, A qR ) — > L 2 {p + An, A qR ) is an isometric isomorphism. The 
inverse of AC is given by the map AC : L 2 {p + An, A qR ) — ► L 2 (p + An, A qR ). 

Proof. Let AC = [AC Mi ^] M ^ e A N be the matrix with elements AC^ ^ = AC(p + p, p + v) = 
p /i ( e i«(p+^)A qR (i/)/((l,l) qR ) 1 /2 an d let A and A be the diagonal matrices with 
the quantities A qR (v) and A qR (/i) iy,p G An) on the diagonal, respectively. (Here 
it is of course again assumed that the columns and rows are ordered by a total 



extension of the partial order (2.3).) We then have that 



AC = A^^KA 1 ' 2 , 



where K — [K ll _ v \ tl _ v ^^ N denotes the matrix given in (4.12). The unitarity of 



the (real orthogonal) matrix K (for parameters in the positivity domain) and the 
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inversion formula K 1 = K imply for the matrix /C that 
(6.13) £*A£ = A and that K.' 1 = t 

(where ft and IC* are the dual and the adjoint (transpose) of JC). Clearly, For- 
mula ( |6.13| ) boils down to a reformulation of the statement in the theorem. (The 
first property in ( 5.13| ) says that the matrix /C determines an isometry between the 



inner product spaces endowed with the positive sesquilinear forms associated to A 
and A, respectively.) □ 

It is instructive to view the g-Racah transform as a Fourier type transformation 
between the grid functions in L 2 (p + A N , A qR ) and L 2 (p + A N , A qR ): 

(6.14a) f(P + »)= (K/X/5 + M) = %==, 

A , 

(6.14b) f{p + V ) = {lC- 1 h{p + v) = -^7==, 

V » ' i /A 

where (■, -)a is taken from ( |6.£| ) and (•, ■) ^ denotes its dual version (and the func- 
tions P M and P v stand for the restrictions to the grids of -P M (e Ma: ) and P^(e lax ), 
respectively). (Recall also that for parameters in the positivity domain ( |6.5a ), 



(3.5b) one has that (1, 1}^ = (1, 1}a as a consequence of Remark 4.1.) 

We will next discuss the behavior of a discretization of the difference operator 
D (|2~4]) with respect to the g-Racah transform. Let L> qR : L 2 (p + A Nl A qR ) — > 



L 2 (p + An, A qR ) be the discrete difference operator of the form 

(6.15) £ qR = £ V +j (p + u)(T j -l)+ V^{p + u){T^-l) 



where 

with 

v(0 = 

AO = 



l<j<n l<j<n 

f+ej£Ajv v— ej-EAjv 



V±j(x) = w(±Xj) J][ v(±Xj + X k )v(±Xj - x k ), 

l<fe<7l 

k¥=3 



•0 



sin(^) ' 

sin f (g, + Q cos f (gt + Q sin | (g e + 1/2 + Q cos f (g d + 1/2 + fl 
sin (M) C0S (M) sin f(l/2 + cos f(l/2 + 



and the action of the operators T ±: is given by 



(T± 1 f)(p + v)=f(p + v±e j ). 

Notice that the conditions v + e j G A at and v — ej G A^v in the summations of 
( 3.15| ) guarantee that the function (Z? qR /)(p + v) for v G Ajy indeed depends only 



on the values of / on the grid points in p + An. Up to (multiplication by) an 
overall constant factor with value ^{tot^hq- 1 ) 1 ' 2 = e^K^-^f+Cso+ ' +ss)/^ 



the operator D qR (6.15) amounts to the restriction of the operator D (2.4) to grid 
functions (cf. Remark B.2) rewritten in trigonometric form (recall ( |6.l[ )). Such 



a restriction is well-defined because the coefficients of the difference operator D 



(p.4|) are regular on the grid points and they vanish (cf. Lemma A. 2) when the 
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q-shift operators T * in (2.4) shift the argument of a function out of the grid 



(which leads us to the above-mentioned restrictions on the sums in the discretized 



operator (6.15)). Let us furthermore introduce the multiplication operator E 



L 2 (p + A N , A qR ) -» L 2 (p + A N , A qR ) defined by 
(6.16) {Ef){p + v) = E v f{p + u) 

with 



l<j'<n 



E v = 2 N ( cosa(pj + Vj) — cos(apj) 



together with the corresponding dual operators D qK : L 2 (p + An, A qR ) — ► L 2 (p + 
Ajv, A qR ) and E : L 2 ( ( o + Ajv, A qR ) -» i 2 (p + A at, A qR ) in which g r -v g r and 
p — > p. (The dual quantities = 2^ 1<<n (cosa(/5j + A 3 ) — cos(aj5j)) coin- 
cide with the eigenvalues £\,a in ( |2-6| ) up to multiplication by the constant factor 
t n-1 (io*i*2*3<T x ) 1/2 = e M [("- 1 )f+(so+ -+g3)/2] relating ^qR ( p7i5[ ) an d £) Q .) 



The following theorem states that the g-Racah transform (i.e the discrete inte- 
gral transformation IC) is the eigenfunction transformation that diagonalizes the 
operators Z? qR and £> qR . 



Theorem 6.5. For parameters in the positivity domain (6.5a), (6.5b) the discrete 
difference operators D qR : L 2 (p + A N , A qR ) — > L 2 (p + A N , A qR ) and D qR : L 2 (p + 
A N , A qR ) — ► L 2 (p+A N , A qR ) are self-adjoint and the map IC : L 2 (p+A N , A qR ) — > 
L 2 (p + Ajv, A qR ) constitutes the corresponding (unitary) eigenfunction transforma- 
tion diagonalizing these operators 

(6.17) K.D qR fC- 1 =E, K,- 1 D qR K, = E. 

Proof. Clearly it is sufficient to prove only one of the diagonalization formulas in 
( 3.17j ) because the other will then automatically follow upon dualization. Since for 
parameters in the positivity domain the (real) functions P\(e la ^ p+!y - ) ), A 6 Ajv form 
an orthogonal basis for the space L 2 (p + A^r, A qR ), it is enough to show that 

(6.18) {KD qK )P x = (EK)Px, for all A <E A N 

(where Px stands for the discretized trigonometric polynomial Px{e la ^ p+ ^)) . Equa- 
tion ( |6.18 ) is immediate from the discretized eigenvalue equation (cf. Remark B.2) 

D^Px = E x Px 

and the observation that (ICPx)(p+p) is nonzero only in the point p+X (i.e. for p = 
A) in view of the orthogonality of the multivariable g-Racah polynomials. (We thus 
have that K.D qR Px = E x ICPx = EKP\.) The self-adjointness of £> qR and Z> qR now 
follows from the fact that the discrete difference operators are unitarily equivalent 
to the real (for parameters in the positivity domain) multiplication operators E 
and E, respectively. □ 

Remark 6.1. If the parameters satisfy the additional constraint 
(6-19) g a - gb - g c - gd = 0, 



then g r = g r , r = 0, . . . ,3 (see (6.4)). Hence, we are then in a self-dual situation (cf. 
Remark O) with both the Hilbert spaces L 2 (p+A N , A qR ) and L 2 (p + A N , A qR ) co- 
inciding and the map JC : L 2 (p + A^v, A qR ) — > L 2 (p + A^, A qR ) being an involution 
(/C 2 =KK. = Id). 
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Remark 6.2. A positivity domain for a self-dual one-parameter subfamily of the 
one- variable g-Racah polynomials with \q\ = 1 similar to the domain considered in 
this section can be found in Section 3C2 of Q together with a discussion of the 
corresponding finite-dimensional discrete integral transform. 



Remark 6.3. The trigonometric polynomials p\(e lax ) are invariant with respect 



to permutations, sign flips (xj 



-Xj) and translations (xj — > Xj + 2n/a) of 



the variables Xi, . . . ,x n . A fundamental domain for R™ modulo the action of the 
discrete symmetry group generated by the permutations, sign flips and translations 
over a period 2n/a is given by the (Weyl) alcove 

(6.20) {x G R"| ir/a > xi > x 2 > ■ ■ ■ > x n > 0}. 

It is therefore natural to consid er th e trigonometric polynomials p\(e zax ), A G A 
as polynomials over the alcove ( 6.20| ). In the trigonometric context the necessity 
to truncate the grid with points p + is, is G A on which the masses of the discrete 
orthogonality measures are concentrated arises as a natural consequence of the pe- 
riodicity of the trigonometric functions (which demands that the support of the 
discrete orthogonality measure be finite). The conditio ns (|6.5a ), (6.5b) in Proposi- 
tion 6.1 arrange things in such a manner that A qB (is) (|6.3j) is positive for v G An 
and zero for v G A\Ajv, and guarantee furthermore that the grid p+An supporting 
the discrete orthogonality measure for the polynomials p\(e lax ), A G An fits in the 



fundamental domain (6.2C) 



Remark 6.4. The second-order operator D (2.4) sits in a commutative algebra gen- 
erated by n-independent commuting analytic difference operators D\ , D 2 , • • • , D n 



of order 2,4,... , 2n, respectively [Dl, D2j. The first of these operators, viz. 
Di, corresponds to the Koornwinder-Macdonald difference operator D of Sec- 
tion |2|. After restriction to functions on the grid p + An the analytic differ- 
ence operators go over in a family of (commuting) discrete difference operators 
Df 1 : L 2 (p + A Nl A qR ) — > L 2 {p + A N , A qR ) (r = 1, . . . , n) given explicitly by (cf. 
Remark B.2 ) 

(6.21) 

Df= U Je , r -\j\{p + v)V e j,jc{p + v)T eJ , r = l,...,n, 

JC{1,... ,n},0<|J|<r 

£j=±l, jeJ; e e j + i/eAjv 

with (T eJ f)(p + u) = f(p + u + e e j), e £ j = J2,eJ e i e i> and 



V £ j,k(x) 



Y[w{e jXj ) Y[ v{e 3 x.j +e r x r )v(e j x j +e r x r + 1) 

3&J 



3<3' 



Y\_ v(ejXj + x k )v(ejXj - x k ), 



keK 



u K , P ( x ) 



+ e v x v )v{-eix l 



H'XV 



LCK, \L\=p l£L 

ei=±l, l£L 



l<l> 



Yl v(eixi + x k )v{em - x k ) 



k£K\L 
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(where v and w are the same as in ( |6.15| )). In the above expressions we have 
used the conventions that empty products are equal to one, and that Uk, p = 1 for 
p = 0. Notice that the coefficient functions V e j. k{x) and Uk, p {x) are regular for 
x G p + Ajv and that the condition e e j-\-v G Ajv in the summation again guarantees 
that (D qR /)(p + v) only depends on the val ues o f / on the grid points in p + An 
when v lies in A^v- (Hence the operator D qR (6.21) is well-defined as an operator in 
L 2 (p+An, A qR ).) After introducing also the corresponding multiplication operators 
E r : L 2 {p + A N , A qR ) — ► L 2 {p + A N , A qR ) given by 



{E r f){p + v) = E r>v f(p + v) 

(r = 1, . . . , n) with 

E TtV = 2 r j2 (- i ) r_iji (ii cosa ^ + ^) 



JC{1,... ,»} 

0<\J\<r 



x ^2 cos(ap ;i ) • • •cos(ap; r _ |J| ) 



r</i<---<^ T ._iTi<n 



(where the second sum in E 1%v should be read as 1 when \J\ = r) and the associated 
dual difference operators Df 1 : L 2 (p + A Nl A qR ) — ► L 2 (p + A Nl A qR ) and dual 
multiplication operators E r : L 2 (p + Ajv, A qR ) — > L 2 ( p + An, A qR ), we are in the 
position to formulate a generalization of Theorem p.5| pertaining to these higher- 
order discrete difference operators. 



Theorem 6.6. For parameters in the positivity domain ( 6.5a ) ; ( |6.5b| ) the com 



muting discrete difference operators Df- : L 2 (p + A N , A qR ) — ► L 2 (p + A N ,A qR ) 
and Df : L 2 (p + A N ,A qR ) — > L 2 (p + A N , A qR ) are self-adjoint and the map 
K, : L 2 (p + Ajv, A qR ) — ► L 2 (p + A^r, A qR ) constitutes the (unitary) joint eigenfunc- 
tion transformation simultaneously diagonalizing these operators 



KDfK,- 1 = E r , 



Kr^DfK, = E r 



(r=l, 



,n). 



The proof of Theorem x6 runs along the same lines as t hat of Theorem 6J3 and 
hinges on the discrete difference equations in Remark B.2 and the real-valuedness 
of the multiplication operators E r and E r for parameters in the positivity domain 
(3.5a), (pT5t ) . For r = 1 Theorem 3.6 reduces to Theorem 3.5. 



Appendix A. Proof for the symmetry of D 



In this appendix we prove Proposition p.3| 



which stated that the g-difference 
•)^ R Km, for 



operator D (2.4) is symmetric with respect to the bilinear form y,-/ N 
parameters satisfying the truncation condition (3.6). This proposition was a key 
ingredient in our orthogonality proof for the multivariable q-Racah polynomials 
(i.e. the multivariable Askey- Wilson polynomials with parameters subject to the 
truncation condition) with respect to the bilinear form (•, ( 3.17 ) (Theorem 3.2). 
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Let us first recall the explicit form of the functions V±j (z) that determine the 
coefficients of the difference operator D (2.4): 

(A.i) y £J (z) = 



n 

l<fc<r. 



[\-tz]z k )(\-tz]z^) 
(l-zU k ){l-zU^) 



e = ±1. 



The symmetry proof for the operator D hinges on two lemmas. The first lemma 
describes a relation between the discrete weight functio n A qR (z/) ( 3.10 ) and the 
coefficients V e j(z) of D evaluated at the grid points rq u (3.7). 



Lemma A.I. Let us assume that v and v + 
(here e is +1 or — I). Then 

(A.2) A« R (v + e ej )V- ej (Tq v+ee * ) = ^ R {u)V ej {rq v ) 



are in the dominant cone A (2.2) 



(where rq" , & qR (v) and V ej (z) are given by (3.7), (3.10) and (A.I)) 



Proof. The relation (|Aj) between A qR = 1/(C* R C^) (|3T0|) and V ±j follows from 
the difference equations 



(A.3) 
and 
(A.4) 



Cf{D) 

Cl K (i> + ej 
Cl R {i>) 



V +] (rq^)f 3 (v) 



for v, v — ej S A 



for D, v + ej e A 



l inkin g the c-functions C^ 1 to the coefficients V±j. The factors fj in (A.3) and 
(A.4) represent certain intermediate products of the form 



/,(„)=* (t Q t lt2 t 3q ) 11 ^ 



l<fc<j 



tTjT^q**-**-!) 



(where we have used the convention that an empty product is equal to one) , which 
cancel each other in the final relation (A.2). The verification of the difference 
equations (A.3), (A.4) is straightforward using the explicit expressions for C+ R , C^ R 
and V±j (in ( |3. 1 1 )- ( |3 . 13 ) and (A.I)) and the elementary shift property (a;q); +1 = 
(1 — aq l )(a; q)i for the g-shifted factorial (I = 0, 1, 2, 3 . . . ). □ 



Lemma A.2. For v G An (3.5) and parameters subject to the truncation condition 
( [3.6| ), one has that 

(A.5) V sj (rq") = if v + ee 3 £ A N (e = +1 or - 1). 



(Here V e j(z) and rq u are taken from (A.I) and (p.7\).) 



Proof. There are two situations that need to be distinguished: either v + eej is in 
A ( |2.2j) but outside the the alcove Ajv (3.5), or v + eej does not even lie in the cone 

The first situation occurs (only) when e = +1 with j = 1 and v\ = N. We 
then have that V+i(rq u ) = because in the numerator the factor (1 — tf,zi) = 



MULTIVARIABLE q-RACAH POLYNOMIALS 



33 



(1 — t n 1 t a tbq N ) becomes identical to zero in view of the truncation 



(l-4ri^ 
condition (|3 

The second situation can occur both when e 
we have that v + eej ^ A iff j > 1 and Vj—i : 
in the numerator the factor (1 — tZjZj_^) = (1 



identically zero. For e 

or if j < n and Vj+\ — i>j. In the former case V. 
contains a factor (1 — t a z~ l ) = (1 — ' ~ 
latter case one has that V-j(rq u ) 

(1 



= +1 or when e = — 1. For e = +1 
Vj. Then V+j{rq v ) = because 

- tT o T 7- 



\qvj-vj-i) = (1 -tt- 1 ) is 
1 we have that v + eej ^ A iff either j = n and v n = 



tzr 1 z j+1 ) = (1 - tTf 1 T j+1 q- v *+"i+i) = (1 _ = o. 



i{Tq v ) = because the numerator 
™) = (1 — t a t~ 1 ) = 0, whereas in the 
because in the numerator one has a factor 

□ 



taT n ^q 



After these preliminaries we are finally set to prove the symmetry relation 
(A.6) (£>/, g)f = (/, Dg)f (/, g € 



for parameters subject to the truncation condition ( |3.6| ). Evidently this amounts 
to showing that (cf. the definition of (•, in (3.17)) 



(A.7) 



J2 V EJ (Tqnf(Tq» + ^)g(Tq»)A^(v) 



l<j<ra,e=±l 



E ^(T/)/(T/)Kr/- ee 0A qR (i>). 



l<j<n,E=±l 



(At both sides of this equation the coefficients V±j are well-defined as rational 
expressions in the parameters subject to the truncati on c ondition (3.6), i.e., no 
den omi nator becomes identical to zero.) Using Lemma A.2 the sums at both sides 
of ( A.7 ) can be restricted resulting in the equation 



(A.. 



E E V ej (rqnf(rq v+ ^)9(rqn^) 



l<j<n i-EAjv with 

e=±l v+eej eA N 



E E V- ej {Tq*)f(Tq>)g(Tq»-^)^{v). 



l<j<n CeA N with 

e=f 1 i>—eej gAn 



To check the identity (A^) (and thus proving (A.6)), one uses Lemma A.l to infer 
that for given j and e each term on the l.h.s. coincides with a term on the r.h.s., 
where v and v are related by v = v + eej. Phrased in other word s: substituting in 
(the terms corresponding to given j and e at) the r.h.s. of (A. 8) v — v + ee,- and 



invoking of Lemma A.l results in (the co rres ponding terms at) the l.h.s. of ( A. 8 ), 
which completes the proof of Proposition 3.3. 



Remark A.l. Notice that in the restriction of the sums, i.e. in passing from (A.7) 
to (A. 8), the truncation condition (3.6) became essential. Moreover, it was only 
after this restriction of the sums that the terms at both sides of the equation 
could be put into one-to-one correspondence. For generic parameters not satisfying 
any truncation condition the proof breaks down at this step because there will be 
nonzero terms of the form 



V +1 (rqnf(rq u+ei )g(rq v )A^) 



34 



J. F. VAN DIEJEN AND J. V. STOKMAN 



with v\ — N in the l.h.s. of ( |A.7| ) that do not match with terms in the r.h.s. and, 
reversely, will there be nonzero terms of the form 

V +1 (Tq>)f(Tqng(Tq^)A^(v) 

with v\ = N that have no counter part in the l.h.s. Therefore, in general equation 
(A. 7) (and hence the symmetry of D with respect to (■, will no longer hold if 
no truncation condition is assumed. 



Appendix B. Computation of orthonormalization constants 

USING PlERI TYPE FORMULAS 

The key to the computation of the sum (p\,p\)ff is a system of Pieri type 
recurrence relations for the renormalized multivariable Askey- Wilson polynomials 

(B.l) P x (z) = Cf(X)p x (z), AeA 

with Cf-(X) given by ( |4.2| ). Basically, the recurrence relations in question provide 
explicit expansion formulas of the type E r {z)P\(z) = c M P M (z) for the products 
of the basis elements P\{z) with certain VK-invariant polynomials E\(z), . . . , E n (z) 
that form a set of generators for the algebra H w of all W-invariant Laurent poly- 



nomials in the variables z\, . . . , z n . Specifically, we have (see [D2, |D3|] ) 



(B.2) E r {z-T) P x (z) = Ujc,r-\J\{Tq X )V £ .,,jc(Tq X )P x+eEj {z), 



JC{1,... ,n},0<)J|<r 

ej = ±l, jSJ; e eJ +AeA 



r = 1, . . . , n, where 

E r (z;r) = ((- i r' Ji n^ + z 7 1 ) 



jc{i,..,t.} je.J 

0<\J\<r 



X 

r</i<---</ 1 ,_i ji <n 



(with r, t taken from ( |3.8| ), flf.q)) and the expansion coefficients are governed by 
V e j, K {z) = II " II Hz- 3 zl*)v(qz £ /zl*) J] v(z £ /z k )v(z £ /z^), 

je.J j,k€J je.J 

j<k keK 

u K , p {z) = {-if J2 (n^H II H^ l 4 k )HQ-^r i ^ k ) 



LCK, \L\=p IGL l,k£L 

e;=±l, IGL l<k 



with 



k£K\L 



-i/2 i'LzK\ = (U,u.n-i\-v» ^<r<^ tr ® 



t)(C) = t- 1 ' A - £ , W(Q = (t tit 2 t 3 q~ 



1-CJ' ^ V™, , (1 _ C 2 )(1 _ gC 2)- 

Here we have employed the notation 
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and we also used the conventions that empty products are equal to one, that Uk. p = 
1 for p = 0, and that the second sum in E r (z; r) is equal to one when \ J\ = r. 

Even though formula (B.2) holds for generic parameters, it is not an entirely 
trivial matter to perform the reduction to the case of parameters satisfying the 
truncation condition (3.6). The problem is that for parameters subject to the 
truncation condition the c- function C+ R (A) ( |3.11 ) becomes infinite for A £ A \ 
Aat. (Recall that we have a zero in the denominator from the factor (tbTi ',q)x 1 = 
(t n ~ 1 i a it,;q)\ 1 = (t n ~ 1 t a tt ) ;q)\ 1 , which is zero for Ai > N when the parameters 



satisfy the truncation relation t n 1 t a tf ] = q N .) Consequently, the renormalized 



polynomials P\(z) ( JB . l| ) are no longer well-defined for such parameters when A lies 
outside the alcove Ajy (3-5). Given the fact that even if we start in the l.h.s. with 
a polynomial associated to a weight A G Ajv, we may end up in the r.h.s. with 
some polynomials corresponding to weights outside the alcove Aat, it is clear that 
in its present form the Pieri type recurrence formulas (B.2) do not make sense for 
all A € Aat when the parameters satisfy the truncation condition ( [3.6] ). 

Lemma B.l. For A G Aat and A + e e j € A \ An, we have that V e j,J a (Tq x ) is the 
product of a factor of the form (l—t a tbt n ~ 1 q) and an expression that is rational in 
the parameters subject to the truncation condition (3.6) (no denominator becomes 
zero after imposing (3.6)/ 

Proof. For A G Aat one has that A + e £ j G A \ Aat iff Ai = N and the index set 
J C {1, ... , n} contains the number 1 with E\ = +1. Then V e jjc{rq x ) picks up a 
factor (1 - i a i b t n - 1 q N ) from the part w(fig Al ) = w(i a t n - 1 q N ). □ 



Lemma B.l tells us that the factor (1 — i a ibt n ~ 1 q N )~ 1 in P\+ eeJ (z) when A 
e £ j G A\ Aa? (stemming from the denominator of the normalization factor C? R (A 
eej)), is compensated in formula ( |B.2| ) by a corresponding factor (1 — i a it,t n 
in the numerator of V e j t j"{rq x ). Moreover, by combining this observation with 



Proposition 3.7 it is seen that if we restrict the variable z to the grid points rq v , 
v G Aat, then we end up with a recurrence formula of the form in (B.2) in which 



the sum in the r.h.s. gets restricted to the weights of the form X + e e j that lie inside 
the alcove Aat. 

Proposition B.2. For parameters subject to the truncation condition (|3.6| ) and 
A, v G Aat one has that 



E r (Tq»;T)P x (Tq») 



E 



Ujc, r -\J\ {rq X ) VeJ, J< (fq X ) Px+e cJ fa"), 



JC{1.... ,n},0<|Ji<r 

=±1, j£J; e eJ +AeA A 



r = 1, 



After shrinking of the domain of the variables to the grid points in the Pieri 
formulas and implementation of the truncation condition, we are now ready to 
compute the the sums (pa,Pa)am A G Aat by means of the same method that also 
led to the norms of the polynomials in the continuous case [D2 . Replacing the 
products at both sides of the identity 

(B.3) {E r P x , Px+„ r )N = {Px,E r P x+ulr y, 



wR 
>N J 



ei H h e r 



by the corresponding r.h.s. of the Pieri formula in Proposition B.2 and using the 
orthogonality of the polynomials P^ with respect to the bracket {■,-) N X leads us to 
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a relation between (P\,P\)'^ L and {P\+ UJr . , m terms of the coefficients of 

the Pieri formula 

(B.4) t/ {1 ,... )r}){r+1) ..., n} (fg A )(P A+Wr) P A+Wr )* = 

f { -i,... ) - r} , {T . + i,...,„ } (fg A+w ")<A,PA>?f 

(recall Uk.p = 1 for p = 0). In the above manipulations we have assumed that A 
and A + to r (i.e. A augmented by the fundamental weight vector uj r — e\ + • • ■ + e r ) 
lie in the cone Ajy and furthermore that the parameters satisfy the truncation con- 
dition. To solve the recurrence relation ( |B.4| ) for (P\, P\)ff we exploit the following 
connection between the c-functions C"J R and the coefficients V{±i,... ,±r},{r+i,... ,n} 

C qR (X) 

+ ^{i,...,r},{r+i,...,«}(^<Z A ) for A, X + aj r eA N , 



Cf(X + U r ) 
C^{\ + LO r ) 



V{-1,... ,-r} : {r+l,...,n}(Tq +UJr ) for A, A + LO r £ A 



JV, 



which is not difficult to derive with the aid the elementary shift property (a; q)i+i = 
(1 — aq l )(a;q)i for the g-shifted factorial (I = 0, 1,2,3. . .). (Notice also that for 
r = 1 these two formulas amount to the (dual versions) of the formulas ( |A.3| ), (A.4) 
in the proof of Lemma A.l specialized to the case j = 1.) Using these two relations 
we can eliminate the coefficient functions Vi±\ .± r }, {r+i n} from (B.4) entailing 

(B.5) (P Al P x )fA* R (\) = (P A+Wr , P x+Ur )fA« R (\ + u r ) 

(where A qR = l/((7+ R C"!i R ) again denotes the 'Plancherel' measure. Since the (fun- 
damental weight) vectors u>i, ... , u> n positively generate the cone A (2.2) it follows 
that the l.h.s. of ( |B.5| ) does not depend on A e An and so we obtain by comparing 
with the evaluation in A = (so P A = 1) that 

(B.6) < P ^ = fe^' AGA - 

which reads in monic form 

(B.7) (px,Px)n -AA qR (A)(l,l) jV AeAjv 

withA^ qR (A) = Cl R (X)/Cf'(X). 

Remark B.l. The proof for the Pieri type recurrence formulas in |D2| is complete 
only for parameters satisfying a self-duality condition of the type 

(B.8) t a q = t b t c t d 

(where t c and td denote the two parameters complementing t a and tb such that 
{t a ,tb,t Cl t c i} — {to, ti, t2, ^3})- This condition on the parameters implies that 
t r = t r and thus that P\(z) = P\(z). As was pointed out in Section 7.2 of [D2|, 
however, the Pieri type recurrence formulas would immediately follow for general 
parameters without self-duality condition once one would succeed in proving that 
Macdonald's conjectured evaluation formula stating that P A (r) = 1 holds for such 
general parameters (cf. also [ D2 , Theorem 3] for a proof of the evaluation formula 
in the self-dual case with parameters subject to the condition ( |B.8| )). At the 'CRM 
Workshop on algebraic methods and g-special functions' in Montreal, May 1996, 
we learned from Prof. Macdonald that he has managed to produce such proof for 
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the evaluation formula with general parameters using an extension of the Cherednik 



approach towards the Macdonald polynomials. (See [CI, C2, M2] for this approach, 
which is deeply connected with the representation theory of affine Hecke algebras.) 
Therefore, we have formulated all our results here without imposing the self-duality 
condition, even though our own direct proof (i.e. without using the representation 
theory of affine Hecke algebras) at present actually requires at one point, viz. the 



verification of the Pieri type recurrence formulas in [ D2 ] (or equivalently the proof 
of the evaluation formula P\(r) = 1), assuming that this additional condition be 
satisfied. 



Remark B.2. Another consequence of the evaluation formula P x (r) = 1 is (sec [ D2 1 ) 
the duality relation for the renormalized multivariable Askey- Wilson polynomials 
originally conjectured by Macdonald 

(B.9) PvlW) = P v (tq»), /i, is E A, 

where P v {z) = C^ R (i l )p l/ (z) denotes the renormalized multivariable Askey- Wilson 
polynomial dual to P v {z) — C'f i (i/)p v (z), i.e. with the parameters t r being replaced 
by the dual parameters i r (|4.l|) . Applying the duality relation to the restricted re- 



currence formulas of Proposition B.2 leads us (up to dualization) to a system of 
discrete difference equations for the multivariable Askey- Wilson/g-Racah polyno- 
mials with parameters subject to the truncation condition 



£ Uje tr _\j\(Tq v )V e j,jc(T<f)Px(Tq u+e '') = E r (rq x -r) P X {rq v ), 

JC{1,-.. ,»},0<|J|<r 

Sj=±l, jSJ; e eJ +!/6Aiv 

r = l,...,n (where V e j,jc and E/j c ,r-|J| are the dual versions of V s j t jc and 
Ujc r —\j\ with the parameters t r replacing t r ). This system of discrete difference 
equations is the restriction to the grid rq v , v G A^ of the system of analytic differ- 
ence equations for the multivariable Askey- Wilson polynomials introduced in |Dl| 



(see also [D2, |D3| ). For r — 1 the discrete difference equation in question boils down 
(after multiplication by a constant factor t" _1 (io ■ • • Hq^ 1 ) 1 / 2 ) to the restriction to 
the grid points of the second order g-difference equation in Section]^ (Eq. Q2.7b )) 
with parameters subject to the truncation condition (p 



£ V ej (r<f) (px(rq v+ ^) -px{Tq v )) - E x , x p x {rq v ), 



l<j<n 

e— ±1; f+eej^Kj^ 



where V s j and E x x are taken from (2.4) and (|2.6| ) 
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